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WG5: Angles from B Decays with Charm 



Charge: 

The main goal is the measurement of the angles y and yS from B decays involving D or char- 
monium mesons. The limitations and ways to overcome them will be discussed, as well as new 
approaches and high statistics projections. 



Introduction 



This report is a summary of the activities of Working Group 5 of the CKM2005 Workshop 
on the Unitarity Triangle (UT), which took place in San Diego, California, from 15-18 March 
2005. The main goal of this Working Group was the determination of the UT angles fi/cpi 
and 7/03 from B meson decays into final states with charm. These determinations play a key 
role for the experimental testing of the Kobayashi-Maskawa (KM) mechanism of CP violation, 
which has already reached an impressive level thanks to the great eff'orts of the BaBar and Belle 
collaborations at the SLAC and KEK B factories, respectively. 

The outline of this report is as follows: after this introduction, a table is provided, listing 
the talks at the San Diego meeting. An executive summary provides a brief overview of the 
activities of the Working Group. This summary describes the current status, both experimental 
and theoretical, and includes an outlook to the future reach of the measurements. More detailed 
information and references can be found in the the collated contributions at the end of the 
document. 




2 



Talks given in Working Group 5 at the San Diego meeting. 



sin(2yS) in Charmonium Modes and Discussion of Experimental Analysis Techniques 



T. Mannel 
T. Browder 
D. Lange 
W.T. Ford 
A. Kusaka 
K. Sumisawa 



M. Verderi 
K. Miyabayashi 
R. Jesik 
T. Gershon 
M. Bruinsma 
B. lyutin 
J. Albert 



S. Ganzhur 
F. Ronga 
C. Voena 
M. Baak 

F. Polci 
V. Sordini 
Y.Y. Keum 

G. Wilkinson 

C. Ferretti 

D. Atwood 
K. Trabelsi 
M. Rama 

Y. Grossman 
A. Petrov 
D. Asner 
R. Fleischer 



J. Zupan 
K. Abe 
M.H. Schune 
N. Neri 
D. Asner 
Q. Zeng 



Theory of the golden mode 

Review of charmonium kaon analysis in Belle 

Review of charmonium kaon analysis in BaBar 

Experimental approaches in BaBar CP measurements 

Vertexing technique in Belle 

Background suppression in Belle 



Measurements related to p 
cos(2yS) with Bd J/iJrK* Experimental Review 
CP inb ^ ccd review 
Bs Jl^<t> status at Tevatron 
Feasibility of jS measurement with Bd Dn^ 
B —> double charm status in BaBar 
Bs —> double charm status at CDF 
7 from B D^*^D^*^ 



Measurements of sin(2jS + y) 
BaBar status and prospects for CP asymmetry measurements 
Belle status and prospects for CP asymmetry measurements 
Extraction of2/3 + y from B D^*'>jt 
B D*p feasibility study 
2yS + 7 from B D'^Ksn' 
2/3 + y from B D^K'-^^^^ 
7 with PQCD 
B, D,K etc. at LHCb 
Ab Jlij/A^ feasibility study at ATLAS 

7 from B DK decay modes 
Combined strategies for y 
Belle status and prospects 
BaBar status and prospects 
Methods to extract y with B DK modes 
Charm input for jS and y measurement, theory 
Charm input for ^ and y measurement, experiment 
7 from Be decays 



7 from B DK with D Dalitz analysis 
Theory introduction 
Belle status and prospects 
BaBar status and prospects 
Statistical treatment in BaBar DK Dalitz analysis 
Experimental Input from CLEO-c 
D Ti^TCn^ feasibility study 



A. Soni 



WG5 summary talk 



Links to these talks are available at http://ckm2005. ucsd.edu/WG/WG5/allWG5 ^essions.php 



3 



Executive Summary 



The leading decay of hadrons containing the b quark is via the Z? ^ c transition. Therefore, 
B decays to final states containing either charm or charmonia are abundant, and provide fertile 
ground for investigations of the phenomenology of the B system, including properties of the 
CKM matrix. For example, the flavour-specific semileptonic decays b civ are used to mea- 
sure B-B mixing properties (in both B^ and 5, systems) such as the mass difference Am and 
the parameter of CP violation in mixing, \q/p\. Furthermore, b ^ c transitions to final states 
consisting of two vector mesons can be used to measure polarization, search for triple-product 
correlations, and, most notably for the decay B^ J/il^(f>, to measure the lifetime difference 
AFv/r,. In addition, the b ^ c transition is the principal tool to investigate the properties of 
heavy baryons and hyperons. While the activity of this working group has touched on these 
areas, the charge restricts the scope to the - nonetheless wide - subject area of measurements 
of angles of the Unitarity Triangle. 

Table 1 lists some of the modes of interest. These modes include measurements of weak 
phases using both mixing-induced and direct CP violation. At the top of the list is the so-called 
"golden mode," B^ J/tf^Kg . This mode allows a theoretically clean measurement of sin(2y0). 
Although there is a penguin contribution, this amplitude has, dominantly, the same weak phase 
as the leading tree diagram, and therefore does not cause significant pollution of the result. 
There has been some recent activity to try to quantify the level of theoretical uncertainty, due 
to subleading terms in both mixing and decay amplitudes, with the outcome that it is below the 
level of the current experimental systematic errors - about 0.01 on sin(2;S) - and likely to be 
much smaller. 



Table 1: Quark level transitions of interest. For each, an alternative transition which gives 
the same final state quarks is given, and whether this alternative is penguin (P), or tree (T) is 
noted. The CKM suppression, including the phase, of the alternative transition is noted, where 

for penguin amplitudes short-distance dominance is assumed (only the top quark in the loop 
is accounted for). Where the alternative is a tree diagram, interference can only occur if the 
ultimately final state can be produced by both intermediate sets of quarks, which contain cu and 
CM, respectively. 

Transition Alternative CKM Suppression Typical decays 

b —> CCS b —> sec (P) 1 Bd —> JI^Ks,Bs —> Jj^cp, . . . 

ccd dcc(F) Rte'^ 5^ ^ 7/i/r;r°, 5^ ^ DW+D^-, . . . 

cud ucd(T) A^RuC^ B^ ^ Dn^ , Ba ^ D^*^^n^ , B, ^ DK^*\ . . . 

cus ucs{T) R^e'y B^^^ ^ DK^*\ B, ^ D(t>, B, ^ Df^K^ . . . 



On the experimental side, this mode is equally golden, with large product branching frac- 
tions and a very clean signal when reconstructed via J/(^ ^ l*l~, Ks n*n~ (I = e,/i). The 
precision of the current results from BaBar and Belle, consequences of the excellent perfor- 
mances of the PEP-II and KEKB accelerators, is such that the systematic errors are no longer 
very small compared to the statistical errors. Detailed understanding of the uncertainties due 
to vertexing, resolution, flavour tagging, and so on, will thus continue become more and more 
important. To reduce these errors, there is particular benefit from the continued operation of 
both B factories, since the two collaborations use somewhat different techniques, and much is 
being learned from discussions and comparisons. The precision of the - already rather mature - 
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sin(2/?) measurements should continue to improve in line with the accumulated luminosity over 
the next few years. In addition, comparisons of sin(2/3) measured with different final states of 
the b CCS transition (eg. J/ij/Ks, rj^Ks , XcoKs) can test some new physics models. 

Measurements of sin(2;0) provide constraints corresponding to two bands in the (p, fj) plane. 
To remove this ambiguity, cos(2yS) should be measured. In general, where there is interference 
between CP-even and CP-odd final states there will be sensitivity to the cosine of the weak 
phase. Several methods based on this concept exist in the literature. For example, in the decay 
Bd JI^K* , where K* Kgn^, diff'erent partial waves contributing to the vector- vector 
final state have different CP properties, and their interference allows measurement of cos(2/3). 
In order to perform this measurement, the relative magnitudes and phases of the contributing 
partial waves need to be determined, which can be achieved with angular analysis. Furthermore, 
to determine the sign of cos(2;S), an ambiguity in the solutions for the strong phase has to be 
resolved. This can either be taken from theory, under certain assumptions, or (as performed 
by BaBar) can be extracted from data by exploiting the interference between K*{S92) Kn 
and the contribution from the broad K-n S-wave in the same invariant mass region of 5 — > 
Jjil/Kn. Results to date prefer the Standard Model solution cos(2yS) > 0, albeit with rather 
large uncertainities. Recently, a new method has been proposed, which utilizes the interference 
pattern in Z) — > KsTT'^tt' decay following DvP (and similar decays). [This method was 

first discussed in WG5 at CKM2005.] In this case the decay model is fixed from studies of 
flavour tagged D mesons (from D*- Dn-) decay, and the ambiguity is removed by using 
Breit-Wigner line-shapes. This method also allows to test the Standard Model prediction that 
the weak phase measured mb ^ cud transitions should be the same as that mb ^ ccs decays 
{ie. 2/3). Results from this method are consistent with the Standard Model, and cos(2yS) < is 
now ruled out with greater than 95% confidence. 

In the limit of tree dominance, and within the Standard Model, B,i decays via the b ccd 
transition also probe sin(2jS). However, the situation is encumbered by the penguin amplitude, 
which is not dominated by a single weak phase (the contributions with t, c and u quarks in 
the loop appear at the same level of CKM suppression, so the assumption of short-distance 
dominance is not well justified). Studies of decays dominated by theb ^ d penguin amplitude, 
may provide more information about this contribution. If there is found to be a non-negligible 
effect due to the weak phase of the penguin, then direct CP violation may arise, allowing an 
additional probe of the weak phases involved. To date, measurements of time-dependent CP 
violation inb —> ccd transitions have been performed using J/i/rn^, D^D~, D*'^D'^ and D*^D*~, 
all of which are consistent with tree-dominance and with the Standard Model. Future updates 
will reduce the statistically dominated errors on these results, and additional channels may also 
be added. 

The decay Bd D^*^-n'^ can proceed either via the doubly-Cabibbo-suppressed b ucd 
transition, or by Bd-Ba mixing followed by the Cabibbo-favoured b cud transition. The 
interference of these two amplitudes results in sensitivity to sin(2y8 -I- 7), and the size of the in- 
terference effect depends on the ratio of the magnitudes of the two amplitudes (usually denoted 
Rd(')k)- Although the D^*^'^n'^ final states are abundant, the smallness of /?D<*);r makes the CP 
violation effect hard to measure and, since it must be extracted from a large number of events, 
sensitive to systematic errors. Furthermore, while there are two observables for each final state, 
there are also two hadronic parameters (Roi'h ^^d (Jdc.);^' the strong phase difference between the 
decay amplitudes), and therefore it is difficult to cleanly extract the weak phase information, 
although approaches based on, eg., SU(3) symmetry exist. A further complication arises due to 
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multiple ambiguous solutions. Nonetheless, these modes are being actively pursued experimen- 
tally. Both BaBar and Belle have results with D-n^ and D*-n^\ BaBar have also investigated 
D-p"^. Both experiments have used partial reconstruction techniques (in addition to the conven- 
tional full reconstruction) to increase the signal yields in the D*"^?:^ channel. Although rather 
different techniques have been employed, most notably to deal with the troublesome possibility 
of CP violation effects on the flavour tagging side of the T(4S ) — > BB event, the results are 
consistent at the current level of precision, which is starting to probe the region where CP vio- 
lation effects are expected to be found. Future updates are therefore of great interest, although 
it will be difficult to continue reducing the systematic uncertainty. 

Some related modes may also prove useful to measure 2jS + y. B4 ^ D*-p'^ has a vector- 
vector final state, and the interfering amplitudes result in an increased number of observables, 
so that all parameters can, in principle, be extracted from the data. However, this mode is exper- 
imentally challenging, eg. the polarization measurement is sensitive to systematic effects, and 
recent results suggest that the CP violation effect may be even smaller than previously expected. 
Larger effects are expected to be found in modes mediated by the b cus & b ^ ucs tran- 
sitions, such as Bd DK^*^ and B^ D^*^*Ksn'^ (note that for time-dependent CP violation 
effects to arise, the kaon must be reconstructed in a strangeness nonspecific final state). These 
techniques are currently limited by statistics and, in the latter case, due to lack of knowledge of 
the resonant structure of the three-body decay. Another interesting mode is B^ D''*^~K^. In 
this case the B^ mixing phase 4>s replaces 2/3, so the time-dependence probes 0^ -1- 7. Since 
is expected to be reasonably large, there will be sufficient observables to extract all parameters 
from the data. The problem of multiple solutions remains, however, this can be addressed firstly 
using additional observables which arise if AFv/F^ is not small, and secondly using U-spin sym- 
metry, which translates s <r^ d, and thus relates 5, D^*^~K^ to B^ D^*^-7t^ . This symmetry 
can be applied in various different ways, and mechanisms exist to quantify the uncertainty due 
to its breaking. This approach looks very promising for the LHCb experiment. 

The decays 5„ j /)(*)^(*) provide the cleanest method to determine 7. The method em- 
ploys the interference between b cus and b —> ucs when the final state is accessible to both 

and mesons. The theoretical uncertainty is completely negligible, and effects due to mix- 
ing and CP violation in the neutral D sector can be taken into account if they are discovered. 
As before, there are important hadronic parameters: the ratio of the magnitudes of the two am- 
plitudes and the strong phase difference between them (these are usually denoted rg and Sb), 
that can be extracted from the data. It is important to emphasise again that this method applies 
to any final state which is accessible to both D° and D^, and one important way to progress is 
to add more decay modes. The theoretical framework exists to do so, both for exclusive and 
inclusive final states. To this end, it is crucial to develop experimental techniques to optimally 
combine the results from the different modes and obtain the best sensitivity on 7. 

The B decay modes which have been exploited to date are B^ DK-, D*K- and DK*"^, 
and in each case the D decay modes to CP eigenstates (principally K^K~ for CP-even and 
Ksn^ for CP-odd), doubly-Cabibbo-suppressed final states (^^tt*), and multibody final states 
{Ksn^n~) have been used. Although fairly large samples of CP eigenstate decays have now 
been accumulated, the statistics are still not sufficient to probe the small CP violation effect 
expected (rg ~ 0(0. 1)). Larger effects are expected using the doubly-Cabibbo-suppressed final 
state; however no significant signals in 5* [K^n-]oK-, and related modes, have yet been 
observed. 
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One of the major developments of the B factories over the past few years has been the use of 
the muhibody decay D Ksn'^n' . The rich interference pattern across the Dalitz plot results 
in regions which are highly sensitive to 7; in addition this mode is reasonably clean experimen- 
tally, due to its large product branching fraction and clean signal. However, in order to perform 
an unbinned fit (necessary to extract the maximum possible information from the data) it is nec- 
essary to make an assumption about the strong phase variation across the Dalitz plot (achieved 
using a sum of resonant amplitudes, parametrized using the Breit-Wigner formalism, and a non- 
resonant term), which results in model uncertainty, currently estimated to be ~ 10° on y. To 
reduce this, additional studies of the Dalitz plot structure are necessary, and the results using D 
mesons coherently produced in —> DD at charm-tau factories play a particularly crucial 

role. When one D meson is tagged as a CP eigenstate and the other reconstructed as Ksn'^n', 
the Dalitz plot density differs from that in the flavour tagged samples, and the differences are 
sensitive to the strong phase at each point. Thus, the D decay model may be verified by fitting 
the CP tagged sample, and if this proves impossible, a model-independent technique (dividing 
the Dalitz plot into bins) can be employed. Additional methods to understand the model uncer- 
tainty (eg. using different multibody D decays) also exist. The understanding of these issues 
has been advanced through informal discussions between members of various collaborations, 
and also theorists; it is expected that this will continue to be a fruitful source of progress. 

By summer 2006, it is to be expected that the available samples at the B factories will 
be approximately double those from CKM2005, and hopefully the combined data sample will 
exceed 1 ab~\ This will allow significant progress in many of the channels discussed above. The 
precision of the sin(2jS) result will continue to improve, and it should be possible to definitively 
rule out cos(2jS) < 0. Significant CP violation effects may be seen in some b ccd channels, 
and also in B^ D^*^-n^. Further progress on the extraction of 7 from 5* £)(*)^(*)± 
expected. Since the precision of the 7 measurement depends on the rg values, it is not possible 
to say with certainty how the error will evolve, but there is some optimism that the CP violation 
effects will cross some threshold of significance. Another exciting prospect is that first results 
from B^ £)(*)^(*) may become available. In additional to other possible new results from 
BaBar and Belle, the output of CLEO-c will be watched closely by those interested in 7. 

Shortly after, attention will increasingly turn to the start up of the LHC. It is clear that de- 
spite (and, perhaps, partially because of) the best efforts of the B factories and the Tevatron 
experiments, there will be potential for LHCb to make a big impact, in particular in the mea- 
surement of Bs mixing parameters and of 7 (in addition to channels outside the scope of this 
working group). As methods exist to measure both /3 and 7 to within very small theoretical 
uncertainty, the next generation of B physics experiments has a promising future. Furthermore, 
these theoretically clean modes act as Standard Model reference points for other channels (eg. 
charmless B decays), with which LHCb and a Super B Factory may search for new physics. 
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1 Time-Dependent CP Violation Measurements in {cc)K^ 

The measurement of the mixing-induced CP asymmetry in the so-called gold-plated mode 
5° J/^ Ks is becoming a precision measurement. This served as a motivation to revisit 
the analysis of this mode within the standard model [1]. 

It is well known [2] that the time dependent CP asymmetry 

a^JJ it) = Cji^ Ks cos(Am t) - S Ks sin(Am t) ( 1 ) 

is very cleanly related to the CKM angle [3, since 

Sjm'Ks ~ sin(2;e), O/vF/f, « 0, (2) 

where deviations from these relations are expected at the level of a few percent. 

The reason for this is summarized as follows. The box diagram shown in fig. 1 leads to a 

A5 = 2 interaction transforming a BP into a B . Thus the Eigenstates of the Hamiltonian are (to 
leading order) the combinations 

p\B')±q\B) q 
|5l/h) = — , with - ^ (3) 

n — " 

where Mn is the matrix element of the box diagram contribution between a B and a B state. 
Due to the CKM factors and the large top-quark mass M12 is completely dominated by the top 
quark. Thus M12 is directly proportional to the CKM phase 

M\2 oc {VtdV;f,f oc and thus - = g-'^^ (4) 

P 

A CP asymmetry is induced by contributions to a decay which have both different weak 
and strong phases. Furthermore, the time evolution also induces a phase difference between the 
two eigenstates 5l and 5h which is exp(?Am t) and which has the same effect as a strong phase 
difference and thus also results in a phase difference. 

The CP asymmetry parameters given in (1) are given in terms of the mixing parameters, the 
decay amplitude A{B^ J/^Ks ) and its CP conjugate as 

r H c 2ImU] 

"YT^ and Sj,^Ks-Y^ 



with 



' = [-pjAiBO^J/^Ks)- 



n.c.t. h d W 

— . « fWVW 



wy yw u,c,t 



AAAAA> 



u, c, t 



I u,c,t d W 

Figure 1: Box diagrams mediating AB = +2 transitions. 
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Analyzing the possible contributions to the decay J/i//Ks one finds that the decay ampli- 

tude contains a tree contribution proportional to V*^K.s which is real in the standard parametriza- 
tion. Furthermore, the leading penguin contribution is also porportional to the same CKM factor 
as the tree, while the remaining penguin contribution (due to the up quark) is suppressed by two 
factors of the Wolfenstein parameter. Thus we have 



leading to (2). 



1.1 Estimating Corrections to the "golden relations" 

Contribution from T. Mannel 

The corrections arise from two sources: (i) corrections to the mixing and (ii) corrections to 
the decay. The corrections to the mixing amplitude due to contributions which do not carry the 
phase exp(2//8) have to be suppressed relative to the leading contribution by a factor m^/m^ ~ 
10""^ due to the GIM mechanism. More precise statements can be made by computing the box 
diagrams of fig. 1 [3, 4, 5]. Since vastly different mass scales are involved, one may make use 
of eff'ective field theory methods in which case one may compute M12 in terms of the mixing of 
the time-ordered product of two A5 = 1 interactions into local A5 = 2 operators of dimension 
8. The dim-8 operator relevant for the contributions with a phase different from ip is given by 
the mixing 



(8) 



which yields contributions with the CKM factor (V*i,Vaf and (V*i,Va)(V*i^Vtd). 

The contributions proportional to shift the imaginary part of the ratio M12/IM12I, which 
will contribute to a deviation from the simple relation (2). One obtains 



Ahn 



\Mn\ 



-(2.08 + 1.23)- 10"^ (9) 



where the uncertainties are due to the input parameters, which are the CKM parameters and the 
masses. 

The second contribution to the deviation from (2) is much harder to estimate. It originates 
from the up quark penguin contribution, which induces a different weak phase into the decay. 
The ratio of matrix elements 

{J/^Ks\H(b ^ uus)\f) WubWJ 
{JI^>Ks\H{b ccs)\B ) \^cb\\ycs\ 

where H{b qqs) denotes the contributions to the effective Hamiltonian due to the b qqs 
transition, is small due to the ratio of CKM factors. The ratio of decay amplitudes appearing in 
the parameter X of (6) can be written in terms of r as 

r = ^ -(\ -lirsmy) . (11) 
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The ratio r is estimated by using an approach based on naive factorization. The up-quark 
loop is evaluated perturbatively and yields [6] 



{JI'VKsWib ^uus)\Bq)= (12) 
-^{JI^Ks ||^ {-sb)y_j, {cc)y + ^ CsT''b)y_^ (cT"c)y^\B,) - lu j + in^ 

where ~ nf-ji^ is an average momentum flowing through the up quark loop and fi ~ mb is the 
typical scale of the process. 

The two matrix elements appearing in (12) are estimated by noting that the tree level efi'ec- 
tive Hamiltonian contains both the color octet and the color singlet matrix elements. It is well 
known that in naive factorization the color octet contribution vanishes and that the prediction of 
naive factorization turns out to be too small. Hence we estimate the color octet matrix element 
by ascribing the difference of the predicted and the the observed rate to this matrix element. In 
this way we get 

\{J/'i'Ks\(sb)(cc)\B'}uct\ = (3.96 + 0.36) • 10'^ MeV' (13) 
\{J/WKs\(sT"b)(cT''c)\f}\ = (1.97 + 0.64)- 10^ MeV^ (14) 

from which we obtain 

Re [r] = (-3.62 + 1.55) • 10"^ , Im[r] = (-4.48 + 1.92) • 10"^ (15) 

where the imaginary part is assumed to be entirely due to the imaginary part of the loop calcu- 
lated in (12). 

The eff'ect on the CP asymmetry is given as 

acp{t) = - (sin(2yS) + AS j/^Ks) sin(Am - ^ cos(Am t) (16) 
and we obtain for the parameter 6 

5 = - (1.02 + 0.75) -10"^ (17) 
The change AS j/^pks in the mixing induced CP asymmetry is given as 



j/xpKs = 2 sin 7 Re [r] cos(2;S) - Aim 
which yields numerically (for siny = 0.86 + 0.12) 



\M 



121 



(18) 



ASj/^Ks = (-2-16 + 2.23) • 10"^ (19) 

In summary it turns out that this naive estimate yields a very small correction to the relations 
(2) from the standard model. However, this statement depends crucially on the way the up-quark 
penguin contribution is estimated, while the calculation of the AB = 2 mixing is quite reliable. 

On the other hand, even if the up-quark penguin contribution is grossly underestimated, the 
standard model corrections to the simple relations (2) are too small to become visible in the 
current 5-factory experiments. 
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1.2 A Method to Measure 0i Using D^*^h^ With Multibody D Decay 

Contribution from T. Gershon 

The value of sin(2^i), where <p\ is one of the angles of the Unitarity Triangle [7] is now 
measured with high precision: sin(20i) = 0.731 + 0.056 [8]. However, this measurement 
contains an intrinsic ambiguity: 20i < — > n - Icpi. Various methods to resolve this ambiguity 
have been introduced [9], but they require very large amounts of data (some impressive first 
results notwithstanding [10]). 

A new technique based on the analysis of Dh^. followed by the multibody decay of 
the neutral D meson, has recently been suggested [11]. Here we use to denote a light neutral 
meson, such as ;7r°, 77, p", co. The modes 5" DcphP, utilizing the same B decay but requiring 
the D meson to be reconstructed via CP eigenstates, have previously been proposed as "gold- 
plated" modes to search for new physics effects [12]. Such effects may result in deviations from 
the Standard Model prediction that CP violation effects 'vab ^ cud transitions should be very 
similar to those observed mb ^ ccs transitions, such as 5° J/iffKs . Detailed considerations 
have shown that the contributions from b ucd amplitudes, which are suppressed by a factor of 
approximately 0.02 [13], can be taken into account. Consequently, within the Standard Model, 
studies of Dcph^ can give a measurement of sin(2<;^i) that is more theoretically clean than 
that from 5° —> J/i/fKs [14, 15]. However, these measurements still suffer from the ambiguity 
mentioned above. 

In the case that the neutral D meson produced in 5" DU^ is reconstructed in a multibody 
decay mode, with known decay model, the interference between the contributing amplitudes 
allows direct sensitivity to the phases. Thus l(f)\, rather than sin(2<^i) is extracted, and the 
ambiguity l(p\ < — » n - 2<pi can be resolved. This method is similar to that used to extract 03, 
using 5* — > DK"^ followed by multibody D decay [16, 17]. 

There are a large number of different final states to which this method can be applied. In 
addition to the possibilities for and the various different multibody D decays which can 
be used, the method can also be applied to B^ —> D*h^. In this case, the usual care must be 
taken to distinguish between the decays D* Dn^ and D* Dy [60]. Also, if is not a 
spinless particle, angular analysis [18] will be required to resolve the contributing amplitudes 
to ^ D*h^. 

A detailed description of the method and the results of a feasibility study are given in [1 1]. 
In summer 2005, the Belle Collaboration has presented the first results on 0i using this tech- 
nique [19]. 

1.3 Measurement of cos 2)8/20i with fi^ ^ JI^K* 

Contribution from M. Verderi 

The measurement of the sign of cos 2/3 provides a direct test of the Standard Model in which 
cosip > 0. The interference of CP-odd (A^} and CP-even (Ao,A^} amplitudes in the decay 

B'^ J/ifrK*'^; K*^^ K^^-nP generates a cos 2/3 contribution in the time- and angle- dependent 
distribution of the decay. This contribution appears as a product of cos 2/3 and cosine of strong 
phases differences. The strong phases themselves are measured up to a two-fold ambiguity 
leading to a sign ambiguity on the cosine of the strong phases, and hence a sign ambiguity on 
cos2;S. 



11 



BABAR resolves the strong phases ambiguity pointing out the Kn 5 -wave, beyond the 
dominant K*{%92) P-wave, in the B Jjif/Kn decay and exploiting the known behaviour of the 
Kn S-P relative phase with Kn mass in the A'* (892) region [20]: among the two strong phases 
solutions, only one leads to the physical behaviour of the Kn S-P relative phase with Kn mass. 

The cos 2p parameter is extracted by a time- and angle-dependent analysis of the 5° —> 
Jlil/K*^;K*^ K^n^ decay. The BELLE and BABAR measurements are shown in Table 2. 
From MonteCarlo studies, BABAR estimates that the non-standard cos 2fi = - Vl -0.7312 = 
-0.67 solution is excluded at 86% CL. 

According to BABAR, the prospect for higher luminosities is that cr(cos 2/3) empirically 
scales as cr(cos2/?) = 7 • (L/fb"')"^^^ for luminosities high enough (> lOOfb"^). The cos2;S 
measurement will be statistically limited for the coming years. 

Table 2: BELLE and BABAR measurements of sin2;8 and cos 2fi with 5° J/i/zK*^; K*^ 
K'^tj^, with both parameters free in the fit and with sin 2fi fixed to world average. 





# events 


sin2;0 


cos 2/3 


BELLE [21] 


354 


0.30 + 0.32 + 0.02 
0.731 (fixed) 


0.31+0.91+0.10 
0.56 + 0.86 + 0.11 


BABAR [22] 


104 


0.10 + 0.57 + 0.14 

0.731 (fixed) 


3.32+J5-^^ + 0.27 

2.72!;;fJ ± 0.27 



2 CP Violation in ^ -> ccd Transitions 

Time-dependent CP violation in B decays governed hy b ^ ccd transitions have been measured 
in B^ -> Jjil/n^, D**D*-, D^D*' and D'D*^ by both Belle and BaBar collaborations. Since 
the tree diagrams of b ^ ccd transitions contain Vcb and Vcd which have no KM complex 
phase, CP violation arises from the complex phase in mixing {i.e. 0i) if neither penguin 
nor other contributions are substantial. These decays are therefore sensitive to the unitarity 
angle yS and complement the precise measurement of sin(2yS) through B^ {cc)K'^^*^ decays. 
The presence of a gluonic penguin contribution with a difi'erent weak phase is expected to 
change the magnitude of the CP asymmetry by not more than a few percent [23]. However, 
additional contributions from non-SM processes may lead to shifts as large as A/3 « 0.6 in some 
models [12]. Interference between SM penguin and tree amplitudes can additionally provide 
some sensitivity to the angle y = arg [ -V^^V*JV^^V*^ ] [24, 25]. 

The potential b ^ d penguin contribution can be probed by other processes. The charge 
asymmetry between 5"^ — ^ J/tf/n^ and B~ J/^n~ is one possibility, although such an asym- 
metry also depends on the strong phases. Replacing cc by ss, the B ^ ^n decay generated by 
b ssd transition solely takes place via penguin transition, so search for such kind of rare B 
decays will give us some information about magnitude of the penguin diagram. 

2.1 CP Violation in 6 -» ccd Transitions at Belle 

Contribution from K Miyabayashi 
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J/ij/Tr is a CP-even final state, and therefore the Standard Model (SM) predicts (in the limit 
of tree-dominance) S jj^jfi = - sin20i Aji^^^ = 0. Results from Belle and BaBar are shown 
in Fig. 2. In the Belle measurement [26], 91 candidates are found from a data sample containing 
152 M BB events. Backgrounds are estimated separately for B J/if/X and combinatorial. The 
purity is estimated to be 86 + 10%. The BaBar measurement uses 88 M BB [27]; 49 events are 
reconstructed in the signal region of which 37 are estimated to be signal. 




Figure 2: (Top) MbdMes distributions from (left) Belle and (right) Babar. (Bottom) At distribu- 
tions and asymmetries. Belle plots are from [26]; BaBar plots are from [27]. 

The obtained results are (see also Fig. 5) 

Belle -0.72 + 0.42(stat) + 0.09(syst) -0.01 + 0.29(stat) + 0.03(syst) 
BaBar +0.05 + 0.49(stat) + 0.1 6(syst) -0.38 + 0.41(stat) + 0.09(syst) 

5° D*^D* has a final state containing two vector mesons, which in general is a com- 
bination of CP-even and CP-odd. The CP-odd fraction can be determined by angular analysis 
of the decay products. Results from BaBar on this analysis are described elsewhere in this 
document. Based on a dataset corresponding to 152 M BB pairs. Belle reconstructed 194 can- 
didate events with 67 + 7% purity in this decay mode [28] The CP-odd fraction is found to be 
0.19 + 0.08(stat) + O.Ol(syst), and it is taken into account in the subsequent fit to obtain the CP 
violation parameters: 

Sd'^d- = -0.75 + 0.56(stat) + 0.12(syst) (20) 
Ad'^d- = -0.26 + 0.26(stat) + 0.06(syst). (21) 

These results are shown in Fig. 3. 
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Figure 3: Results of the Belle analysis of 5° D**D*- [28]. (Left to right) Distributions of 
Mbc, cos 0tr and raw A? asymmetry separately for (top) poor and (bottom) good quality tags. 

The final states D*-D^ are not CP eigenstates, and hence the phenomenology is more com- 
plicated, though the time-dependent observables are still sensitive to 20i. Results from BaBar 
on this analysis are described elsewhere in this document. Belle has reconstructed this decay 
mode with two different methods, full and partial reconstruction [29], using a data set of 152 M 
BB pairs. The results are summarized in Fig. 5. 

All measurements mentioned above are statistically limited, so finding new b ccd decay 
modes is also important. Both Belle [30] and BaBar [31] collaborations have observed 5'' 
J/il/p^ decays; this final state consists of two vector mesons and therefore angular analysis is 
necessary to extract the CP composition. Furthermore, the contribution of non-resonant 5" 
Jlif/n^n" decays must be well understood. Recently, the doubly charmed mode 5" D^D' 
has been observed by Belle [32]. This is a CP eigenstate, and looks quite promising for time 
dependent CP violation measurements in the near future. 

2.2 Time-Dependent CP Asymmetries in z)(*)± /)(*)+ at Babar 

Contribution from M. Bruinsma 

Using (232 + 3) x 10^ 7(45 ) ^ BB decays recorded by the BaBar detector [33] at the PEP-II 
e^e' collider, Babar has measured CP asymmetries in 5° decays to the CP eigenstates D*^D*' 
and D'^D~ and CP 'pseudo'-CP-eigenstates D*^D^ and D*~D^. The results are summarized in 
Table 3 and are consistent with S = - s'm{2/3) and C = 0, expected in the SM for a tree-level- 
dominated transition with equal rates for 5° D*^D' and 5" D*'D^. 

The vector-vector final state D**D*~ is predominantly CP-even, but contains also CP-odd 
contributions. The CP-odd fraction for this decay is measured in a time-integrated angular 
analysis to be R, = 0.124 + 0.044(stat.) + 0.007(syst.). The S coefficient for D*+D*- in table 
3 is corrected for this CP-odd dilution. The results for D*^D*~ are preliminary while those for 
D*-D+, D*+D- and D+D'are final[34]. 

2.3 Bs double charm status at CDF 

Contribution from B. lyutin 
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Sample 



Purity 



5 



C 



D**D 



396(23) 



0.75(2) 



0.04 ± 0.14(stat.) ± 0.02(syst.) 



-0.65 ± 0.26(stat.)!^-^^(/?r) ± 0.04(syst.) 
£)*+£)- 145(16) 0.49(3) -0.54 ± 0.35(stat.) ± 0.07(syst.) 0.09 ± 0.25(stat.) ± 0.06(syst.) 
D*-D+ 126(16) 0.49(3) -0.29 ± 0.33(stat.) ± 0.07(syst.) 0.17 ± 0.24(stat.) + 0.04(syst.) 
D+D- 54(11) 0.37(6) -0.29 ± 0.63(stat.) ± 0.06(syst.) 0.11 ± 0.35(stat.) ± 0.06(syst.) 



Table 3: Signal yield and purity for each of the 5" /)(*)± /)(*)+ samples, and results for the 
parameters S and C describing the time-dependent CP-asymmetries. The purity is defined as 
the fraction of signal events in the region niss > 5.27 GeV. 

With 243 pb~^ of CDF displaced track trigger [35] data the ratio of branching fractions 
for fully reconstructed hadronic modes is measured Br(^ —> D'^D~)/Br(B^ —> D~3n) where 
^ Knn and D, (pn{K*K, nnn). 

The optimization of selection cuts was done maximizing the significance of the signal Monte 
Carlo over combinatorial background from data. The signal region was blinded. 

Monte Carlo simulation is used to determine shapes describing signal and partially recon- 
structed decays. The shapes were fixed while doing the mass fit (figure 4). Only relative nor- 
malizations were allowed to float. The mass of the signal was also allowed to vary as well as 
parameters describing combinatorial background. Signal yields were obtained integrating over 
a ±2.5cr window, where a means the fitted signal width. 

The analysis efficiency was extracted from signal Monte Carlo. To reduce systematic Monte 
Carlo mass distribution was fitted and integrated it in a 2.5cr window the way it was done for 
data. 

Among the systematic uncertainties an important one is the 5" D 3n sample composition 
uncertainty. It arises due to the poor knowledge of the 3n sub-resonant structure. Cut efficiency 
uncertainty is assigned due to the difference between data and Monte Carlo distributions for 
variables we cut on. A fit uncertainty is due to the fit instability with respect to fixed parameters. 
For D~D^, where —> 3n, there is also a sample composition systematic due to the poor 
knowledge of sub-resonant structure. All of these combined are well below the statistical 
uncertainty. 

The ratios of branching fractions corresponding to the three different modes: 



Br(fi«^ D-37:) 



1.95 + 0.20(5^0 + 0.12(syst)±0A9(BRi) 



Br{B^^D-2,n) 



1.83 + Q22{stat) + Q.n{syst)+QA6{BRi) + QA1{BR2) 



BriB"^D,D^,D,^7Tnn) 
Br{B^^D-3T:) 



2.46 + 0.34(5^0 + QAl{syst)+ 0.62 {BR^) + 0.34 {BR^} 



Branching fraction uncertainties are quoted separately: 



• BRi is due to Br{Ds (pn) 



• BR2 is due to 



Br{Ds^K*K) 
BriDs^4>n) 
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Figure 4: Data fits of 5° ^ D'Stt and 5° D^D-,D^ 
and 5° ^ D^,D~, nnn (bottom). 



(pn (top), 5° ^ D^D-,D^ K*K 
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• BRi. IS due to „ ' — tt 

where branching fractions from PDG (2004) were used. In combining results several assump- 
tions were made. The systematic uncertainty for all three modes was considered fully corre- 
lated, and "BR2" and "BRi," uncertainties were assumed to be non-correlated. The final value 
is: 



5r(5° ^ D^D ) 
Br(B^^ D-3n) 



2.00 + 0.16(A^C) + 0.12(^3;^0 + 0.50(5i?i) 



where only one number is quoted for the non-correlated error (marked "NC"). The current 
world average from the PDG (2004) is: 



Br{B^ D]D ) 
Br{B^ D-3n) 



1.00 + 0.39 



which is consistent with CDF result within 2cr. The error on (pn branching fraction 

(3.6% + 0.9%), which should be greatly reduced in coming years, dominates the measurement 
for now. This analysis prepares the way to measure Br{Bs D s) I Br{B^^ D^D ). 

To summarize this section, time dependent CP violation inb ^ ccd transitions have been 
measured using 5" J/il/n'-^ and doubly charmed modes. The measurements performed so far 
are summarized in Fig. 5. All these measurements are statistically limited. Note that no direct 
CP violation has been observed so far. With 5 x 10^ or 10^ BB, we would expect to be able to 
observe non-zero CP violation, and in order to check the deviation from b ccs modes, much 
more data is needed. Therefore, it is interesting to continue to hunt for new decay modes. In 
order for better understanding of the b ^ d penguin contribution, it is also important to search 
for rare decay modes and CP asymmetries in relevant charged B decay modes. 



BABAR 

0.722 ±0.040 + 0.023 

Belle 

0.728 ±0.056 ±0.023 

Average (incl. CDF and LEP) 

0.726 ±0.037 



J/>|;7i" (BABAR) 
-0.05±0.49±0.16 

J/i|/ji° (Belle) 

0.72±0.42±0.09 

(B^^ftAR) 

0.65±0.26_,„, 

(Belle) 

0.75±0.56±0.12 

D**D" (BABAR) 
0.82±0.75±0.14 
□•"D* (BABAR) 
0.24±0.69±0.12 

(Belle) 
0.55 ±0.39 ±0.1 2 

□•"D* (Belle) 

0.96±0.43±0.12 



BABAR 

0.051 ±0.033 ±0.01 4 

Belle 

-0.007 ±0.041 ±0.033 

Average 

0.031 ±0.029 

J/H/n" (BABAR) 
0.38±0.41 ±0.09 

ilv/n" (Belle) 

0.01 ±0.29 ±0.03 

D*" (BABARr" 
0.04±0.14±0.02 

□•*D'" (Belle) 

0.26±0.26±0.06 

□•*D" (BABAR) 
-0.47 ±0.40 ±0.1 2 
□•"D* (BABAR) 
-0.22 ±0.37 ±0.10 

□•*D" (Belle) 

-0.37 ±0.22 ±0.06 

□•"D* (Belle) 

0.23 ±0.25 ±0.06 



-Tl xS, 



-1 -0.5 

C = (1-W')/(1 



Figure 5: Summary of time-dependent CP violation measurements m b ^ ccd transitions. 
(Left) -rifS y„ where 77/ is the CP eigenvalue of the final state. (Right) C/^^ = -Af^^. 
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2.4 7 from B D^*^D^*^ and B D^D^*^ 

Contribution from J. Albert 

One can combine information from B £)(*)£)(*) and B branching fractions, 

along with CP asymmetry measurements in 5 ^ D^*^D^*\ to obtain information on y [24, 25] 
(see also [36]). An extraction of constraints on y may be obtained by using an SU(3) relation 
between the B D^*^D^*^ and B D^D^*^ decays. In this technique, the breaking of SU(3) 
is parametrized via the ratios of decay constants /^c*) ///)(•), which are quantities measured in 
lattice QCD [37]. The quality of the constraint on 7 using this technique is highly dependent 
on the difference in the values of the CP asymmetries in 5 ^ /)(*)£)(*) from the value of sin 2fi 
as measured in charmonium decays, as this difference is what provides information on the 
additional weak phase. 

3 Measurements of sin(2)S + 7)/sin(20i + ^3) 

The decay modes jB" D''*^^n- have been proposed to measure sin(2jS+7) [38]. In the Standard 
Model the decays B^ D^*^*n~ and B — > D'^*^'^n' proceed through the b Tied and b ^ c 
amplitudes A„ and A^, respectively. The relative weak phase between these two amplitudes is 7. 

When combined with E^B mixing, this yields a weak phase difference of 2jS + 7 between the 

interfering amplitudes. 

The decay rate distribution for B — ^ D''*^~Tt^ is described by an equation similar to 1, where 
the parameters C and S are given by 

C^l--^, 5*^-^sin(2;S + 7±5). 

1 + r"^ 1 + r"^ 

Here 5 is the strong phase difference between A„ and Ac and r = |A„/At |. Since A„ is doubly 
CKM-suppressed with respect to A^., one expects r to be small of order 2%. Due to the small 
value of r, large data samples are required for a statistically significant measurement of S . 

3.1 Belle IMeasurements of sin(20i + ^3) 

Contribution from F. Ronga 

Belle has probed the CP- violating parameter sin(20i +^3) using full reconstruction of 5° 
D^*^n decays [39] and partial reconstruction of B^ —> D*n decays [40]. Both analyses are based 
on a sample of 152 million BB pairs. 

Belle expresses the time-dependent decay probabilities for these modes as 

P(B^ D^*'^^7i^ « J_e-|Afl/rso [1 + cos(AmAO - 5 * sin(AmAO] 
8tbo 

P(50 ^ £,(*)±^T « -^e-'^"/^5''[l ±cos(AmA0-5*sin(AmA0] (22) 
8tbo 

where S - = (-l)^2i?£,(.)^sin(20i + ^3 ± ^ow^r)- ^ is the angular momentum of the final state (1 
for D*n) [15], Roi'^n is the ratio of magnitudes of the suppressed and favoured amplitudes, and 
5jy(*)„ is their strong phase difference. It is assumed that Rjy{*)„ is small [13]. 
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In the full reconstruction analysis, the following decay chains are reconstructed (charge con- 
jugate modes are included): B D*^n ,D*^ D^n^^,D^ K'tt^ , K-n*n^,D^ K'n^n^n- 
and B D'^n'',D* K'n^n^. The event yields are summarized in Table 4. 



Table 4: Number of reconstructed candidates, selected candidates (after tagging and vertexing) 
and purity, extracted from the fit to (A£', Mbc). 

The standard Belle tagging algorithm is used to identify the flavour of the accompanying 
B meson. It returns the flavour and a tagging quality r used to classify events in six bins. The 
standard Belle vertexing algorithm is then used to obtain the proper-time diff"erence A?. A fit to 
the A? distribution of events in the (Afi, Mbc) side-band region is performed to get the parameters 
of the background PDFs. A lifetime fit is the performed on all components together to get the 
resolution parameters; the wrong-tag fractions for the six r bins are also determined from the 
data. 

The S - have to be corrected to take into account possible tag-side interference due to tagging 
on 5" DX decays. Effective corrections {i'^g}'^** are determined for each r bin by a fit to 
fully reconstructed D*£v events, where the reconstructed side asymmetry is known to be zero. 
Finally, a fit is performed to determine 5", with Am and Tgo fixed to the world average, and the 
wrong-tag fractions and {■S'^g}'^^ for each r bin fixed to the values determined previously. We 
obtain: 



The systematic errors come from the uncertainties of parameters that are constrained in the fit 
and uncertainties on the tagging side asymmetry. The result of the fit for the sub-samples having 
the best quality flavour tagging is shown on Fig. 6. 

The partial reconstruction of 5 ^ D*'^nf, D*^ D'^tTs is performed by requiring a fast pion 
TTf and a slow pion Ks, without any requirement on the Z)°. The candidate selection exploits the 
2-body kinematics of the decay using 3 variables: the fast pion CM momentum p* ; the cosine 
of the angle between the fast pion direction and the opposite of the slow pion direction in the 
CM (cos Sfs); the angle between the slow pion direction and the opposite of the B direction in 
the D* rest frame (cos ^hei)- Yields are extracted from a 3D fit to these variables (see Table 5). 
The flavour of the accompanying B meson is identified by a fast lepton, ^tag- The proper time 
At is obtained from the z coordinate of ;r/ and ^tag constrained to the 5-lifetime smeared beam 
profile. 

The resolution function is modeled by a convolution of three Gaussians whose parameters 
are determined by a fit to a J/i/r jx^jT sample selected the same way as the signal sample. 
In order to correct for possible bias due to tiny misalignment in the tracking devices that would 
mimic CP violation, the mean of the Gaussian resolution is allowed to be slightly offset. The 
vertex position is also corrected to account for possible misalignments. A fit to events where the 



Decay mode Candidates Selected Purity 



B^Dn 9711 9351 91% 

B^D*n 8140 7763 96% 



2Rd„ sin(20i -1-^^3-1- Son) 
IRon sin(20i +(p3- don) 
lRD*n Sin(20i -I- 03 + SD*n) 
IRo'n Sin(20i -I- 03 - do'n) 



0.087 + 0.054(stat) + 0.018(syst) 
0.037 + 0.052(stat) + 0.018(syst) 
0.109 + 0.057(stat) + 0.019(syst) 
0.011 + 0.057(stat) + 0.019(syst) 



(23) 
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Mode 


Data 


Signal 


D*p 


Corr. bkg 


Uncorr. bkg 


SF 


2823 


1908 


311 




637 


OF 


10078 


6414 


777 


928 


1836 



Table 5: Fit yield for the signal and the various types of background in same-flavour (SF) and 
opposite-flavour (OF) events. 



two pions have same sign is performed to determine the shape of uncorrelated background A? 
distribution, while a fit to the {p*p cos ^hei) side-bands provides that of the correlated background. 
The D*p component is modeled the same way as the signal, with S- fixed to zero. 

A fit for Am and Tgo is performed to check the fit procedure. A fit to a D*£v sample selected 
similarly to the signal sample is performed to check the bias correction. Several MC samples 
with S- are fitted to check possible fit bias in the extraction of 5 

An unbinned maximum likelihood fit with Am and r^o fixed to the world average, and S-, 
At offsets and wrong-tag fractions floated yields: 



IRo-n sin(20i + 03 + do'^r) = -0.035 + 0.041(stat) + 0.018(syst) 
IRo'n sin(20i + 03 - do'n) = -0.025 + 0.041(stat) + 0.018(syst) 



(24) 



The main systematic errors come from the background fractions, the background shapes, the 
resolution function and the offsets. The result of the fit is shown as asymmetries on Fig. 7. 
Asymmetries are defined as: 



iSF 



lOF 



Increase of the available data and addition of more modes in the full reconstruction, as well 
as tuning of the selection and vertexing on more Monte Carlo and data, will help reduce both 
statistical and systematic errors. A reduction by a factor 0.3 for the former and 0.5 for the latter 
is expected with 1 ab"^ 




-5 
At(ps) 



5 5 
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Figure 6: A? distributions for the Dn events (top) and D*n (bottom) events with the best quality 
flavour tagging. 



20 



>< 






1 .0 


0) 




E 




E 


1 


>. 


w 




< 






0.5 



<~JZ> 

Belle 



-H — m 

-0.5 
-1 



'0.1 -0.05|^ 




-0.2-0.1 5-0.1-0.05 0.05 0.1 0.1 5 0.2 

Az/cm 



-0.2-0.15-0.1-0.05 0.05 0.1 0.15 0.2 

Az/cm 



Figure 7: Results of the partial reconstruction fit shown as asymmetries for (left) SF events and 
(right) OF events. The inset plots magnify the region around Az = 



3.2 Measurements of sin(2>S + y) at BaBar 

3.2.1 Status and Prospects for CP Asymmetry Measurements 

Contribution from S. Ganzhur 

Two different analysis techniques, full reconstruction [41] were used for and partial recon- 
struction [42] the sin(2jS + y) measurement with 5° D^*^"^!:-. 

The full reconstruction technique is used to measure the CP asymmetry in 5° D^*^^n- 
and 5° D*^p- decays. The result with 1 10 million BB pairs is 

Ir^'' %m{2l5 + y) cos S'^" = -0.032 + 0.03 1 + 0.020 

2r^'^cos(2;e + 7)sin5^'^ = -0.059 + 0.055 + 0.033 

Ir^'"" sin(2;e + y) cos S^'"" = -0.049 + 0.03 1 + 0.020 

2r°*"cos(2y8 + y)sin5°*'' = +0.044 + 0.054 + 0.033 

2r^^ sin(2y8 + y) cos c^^'' = -0.005 + 0.044 + 0.021 

2r^^cos(2y8 + y)sin5^^ = -0.147 + 0.074 + 0.035, (25) 

where the first error is statistical and the second is systematic. The systematic error for 5" 
D*'^p- includes the maximum bias of asymmetry parameters due to possible dependence of 
r on the tttt" invariant mass. For the measurement of 2rcos(2/3 + y) sin 5 parameter only the 
lepton-tagged events are used due to a presence of tag-side CP violation effect [43]. 

In the partial reconstruction of a 5" D*^n- candidate, only the hard (high-momentum) 
pion track tt/j from the B decay and the soft (low-momentum) pion track Ks from the decay 
D*~ CPn~ are used. Applying kinematic constraints consistent with the signal decay mode, 
the four-momentum of the non-reconstructed, "missing" D is calculated. Signal events are 
peaked in the irii-aiss distribution at the nominal D° mass. This method eliminates the efficiency 
loss associated with the neutral D meson reconstruction. The CP asymmetry independent on 
the assumption on r* measured with this technique by BABJ^ using 232 million produced BB 
pairs is 

Ir^'" sin(2;8 + y) cos 6^'" = -0.034 + 0.014 + 0.009 

2r°*'^ cos(2y8 + y) sin S'^'" = -0.019 + 0.022 + 0.013, (26) 
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Figure 8: Raw asymmetry for (a) lepton-tagged and (b) kaon-tagged events of B*^ 
decay mode using the method of the partial reconstruction. The curves represent the projections 
of the PDF for the raw asymmetry. 



where the first error is statistical and the second is systematic. This measurement deviates from 
zero by 2.0 standard deviations. Figure 8 shows the raw, time-dependent CP asymmetry 

NBo(At)-N^m 

A(At) = 

NBo(At) + N-^(At) 

In the absence of background and with high statistics, perfect tagging, and perfect At measure- 
ment, A (AO would be a sinusoidal oscillation with amplitude 2r sin(2;S -I- 7) cos 5. 

Recently it was proposed to consider the D'^*^'^a* 2) decays for measurement of sin(2yS-l- 
7) [44]. The decay amplitudes of B mesons to light scalar or tensor mesons such as or 
emitted from a weak current, are significantly suppressed due to the small decay constants fa^^^y 
Thus, the absolute value of the CKM-suppressed and favored amplitudes become comparable 
and the CP asymmetry in such decays is expected to be large. However, the theoretical predic- 
tions of the branching fractions for 5*^ D^*^'^af^^^^ is expected of the order of (1^4)- 10"^ [45]. 
One way to verify the expectations and test a validity of the factorization approach is to measure 
the branching fractions for the more abundant decay modes B^ — > Z)^*^'^ao(2)- Using a sample 
of about 230 millon Y(4S ) —> BB no evidence for these decays were observed. This allowed 
one to set upper limits at 90% C.L. on the branching fractions B(B^ —> D'^a^) < 4.0 x 10"^ and 
5(5° ^ £)>-) < 2.5 X 10"^ 

The decay modes B have been proposed for determination of sin(2;8 + 7) from 

measurement of time-dependent CP asymmetries [46, 47, 14, 15]. In the Standard Model the 
decays of 5° and B mesons into final state D^*^"^^ proceed through the b ^ c and b ^ u 
amplitudes, respectively. Due to relatively large CP asymmetry (r = \A{B 
/)(*)o^o-j| ^ Q 4-) these decay channels look very attractive for such a measurement. Since the 
parameter r can be measured by fitting the C coefficient in time distributions, the measured 
asymmetry can be interpreted in terms of sin(2yS -I- 7) without additional assumptions. However, 
the branching fractions of such decays are relatively small ~ 5 • 10~^ and a large data sample is 
therefore still required. Moreover, the B decay dynamics can lower the expectation for the ratio 
r. The magnitude of this ratio can be probed by measuring the rate for the decays B £)(*)o^*o 
and B — > D^*^^K*^. From the measured branching fractions [48], one obtains r < 0.8 at the 
90% C.L. from a central value of r = 0.4 + 02(stat.) + O.lisyst.) 
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3.2.2 2)8 + 7 from D^K^n- decays. 



Contribution from F. Polci, M.-H. Schune & A. Stocchi 

The use of the 5" D^K^n decays for the extraction of2/3 + y has been proposed in [49]. 
In these decays the sensitivity to the weak phase 2;S+7 comes from the interference between the 
Cabibbo allowed and the Cabibbo suppressed amplitudes leading to the same final state through 
the B^-B^ mixing. Here the main advantage comes from the possibility of performing a Dalitz 
analysis, which allows to reduce the eight fold ambiguity in the determination of2fi + y to only 
a two fold ambiguity [50]. 

To explore the potentiality of this approach, a study of the sensitivity to 2;S + y has been 
performed. The distribution of the invariant masses Mj^o^- and M^+j^ has been parametrized 
with a model where the decay amplitude in each point / of the Dalitz plot is a linear combina- 
tion with complex parameters of resonances described by Breit-Wigner functions. This model 
realistically reproduces the distribution obtained by the BaBar collaboration [51]. 

The time evolution can be written as: 

|^.|2 _ 1 2\m(r ) 

P(50 ^ D^f^TT-) = ^ "' e-^'{l + ( , ' J cos(AmO + ( , ,/ ' ) sm(Amt)} (27) 

2 Ir^P + 1 |ri|2 + 1 

where A^. (A„.) is the magnitude of the Cabibbo allowed (suppressed) amplitude and r, the ratio 
between the Cabibbo suppressed and the Cabibbo allowed amplitudes which varies across the 
Mko^—Md+ko plane. Analogous expressions can be written for —> D*K^n~, B^ —> D'K^n* 
and B'^ D-K^n* decays. 

. The regions showing the highest sensitivity to 2/?+y are the ones with interference between 
5" D**Of^o and the Cabibbo suppressed 5" 0**0^0 between 5<' D^K*' and the 
Cabibbo suppressed 5^ ^**o^o depends on the actual Z)**° component. 

The conclusion of a feasibility study is that with 500 fb~^ the relative error will lie between 
25% and 50%. 

3.2.3 sin(2/8 + 7) constraint from CP asymmetries in fi** -» D^*^7r/p decays 

Contribution from C Voena 

The CP parameters extracted from the time-dependent evolutions of 5" — > D^*^nlp decays 
that have been studied at the 5-factory experiments Babar and Belle can be written as: 

ttj = 2rj sin(2/? -I- y) cos 6j , 
bi = 2r-sm(2/3 + y)cosS'i , 

Cij = 2cosi2fi + y){rjSm6j-r'.sm6'i). (28) 

where i is the tagging category and j is the reconstructed B decay (j=Dn, D*n, Dp). The param- 
eter rj is the ratio of the suppressed over the allowed amplitude contributing to the corresponding 
decay and 6j is a strong phase. The primed parameters are the corresponding quantities related 
to CP violation on the tagging side. We expect a very small asymmetry in these decays, rj is 
expected to be ~ 0.02 (with ^ < rj). In the extraction of sin(2;8 -1- y) we make use of the CP 
parameters free of the unknown tag-side interference: 

Uj = 2rj sin(2yS -I- y) cos 6j , 
ciepj = 2cos(2p + j)rjsm6j . 
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For each mode we have two observables and three unknowns, external inputs are therefore 
needed to extract sin(2/3 + y). SU(3) symmetry is currently used to estimate the rj parameters. 
The relation (for 5" D'^n~) is: 



r.. = .an«e^| ^g:g;;;:j |[38] (29) 

Similar relations have been used for the other two decay modes. Equation 29 has been ob- 
tained with two approximations. In the first approximation, the exchange diagram amplitude 
contributing to the decay D^n~ has been neglected and only the tree-diagram amplitude 
has been considered. No reliable estimate of the exchange term for these decays exists although 
there are experimental hints that it suppressed with respect to the tree term. The second ap- 
proximation involves the use of the ratio of decay constants fo/fo, to take into account SU(3) 
breaking effects and assumes factorization. We attribute a 30% relative error to the theoretical 
assumptions involved in obtaining the value of rj. 
Using the current experimental inputs we obtain: 

ron = 0.020 + 0.003 + 0.006(?/j.), , 

ro., = o.oi5:;;i;;;^±o.oo5(?/i.), 

rop = 0.006 + 0.003 + 0.002(?/z.) . 

It has been suggested by theorists [52] to use the decay D'^n'^ to determine rnj^ (and 

similar for the other two modes) making use only of isospin symmetry (and not of SU(3) sym- 
metry). The measurement of these branching ratios seems however challenging since we expect 
few events even at high luminosity (~ 6 in 500 fb~^). Using the current world averages from 
HFAG for the aj and ciepj parameter and the rj parameters obtained above we determine con- 
straints on sin(2;S -i- y) [ 1 17] using two different statistical approaches: a bayesian [53] approach 
which gives: 

2fi + y = (SS:ty (30) 
and a frequentistic [54] approach which gives: 

I sin(2;e + y)| > 0.49 (0.27) at 68% (90%) CL. 

Figure shows the corresponding bounds obtained in the unitarity (p, ff) plane. 

Projecting the error on 2 ab~', assuming the same measured central values for the Uj and 
ciepj parameters and no improvement on the knowledge of vj gives an error on 2yS -i- y of about 
30". 

3.2.4 sin(2/8 + y) from fi» ^ D»/i:<»/+> decays 

Contribution from V. Sordini, M. Pierini, L. Silvistrini &A. Stocchi 

In terms of the Operator Product Expansion, we can express the amplitudes for the decays 
of the B to neutral Z)-charged K: 

A (B^ D^K^) = T + C 
A (B^ D^K*) = C + A 
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Figure 9: Bound in the unitarity (p, 77) plane from sin(2jS + 7) in the frequentistic (left) and in the 
bayesian approach (right). For the bayesian approach, the blue region corresponds to the 68% 
confidence level and the gray region to the 90% confidence level. 



where T stays for "tree contribution" (T = V^s V*^ £'1 (s, I, c,K,D)) , A for "annihilation 
contribution" (A = VusVl^Ai{s,l,c,K,D)) , C and C are "colour-suppressed contributions" 
( C = y,, y*^ E2 (c, /, s,D,K),C = y„, y;^ E2 (/, c, s, D, K) ) The parameters Ei ,£2 and Ai are 
renormalization scheme and scale independent. We define roo^+ as the amplitudes ratio: 

_ \A (5+ ^ D^K-')\ 
~ |A (5+ ^ D^K-^)\ 

Similarly for the decays of the B to neutral D-neutral K we can write the amplitudes in terms of 
the same parameters : 

A (5+ ^ D^K^) = C 
A (5+ ^ D°K°) = C 



and define the amplitudes ratio: 



\A (50 ^ D^K^)] 



C is known since the BR(B^ D^K^) is measured and writing the two ratios in terms of the 
OPE parameters 



\C + A\ 
\T + C\ 



\C\ 

ro^KO = - 



we get for C the equation: 

|C| = -A cos 6 + ^A2(cos2 6-l) + r]^^^BR{B^ D^K^) 
where the phase space term cancels and 6 is the strong phase . 
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There are two solutions for |C|, these are both acceptable only if: 



cos(5) < and > rlof.,BR(B^ D^K*) 

If we use the actual experimental values: 

BR(B^ '^K^) = (3.7 + 0.6) • 10"^ 
^ D^K^) = (5.0 + 1.4) • 10"^ 
BR(B* D^K^) < 5.0 • 10"^ @ 90% probability 

( the BR{B^ D"^^'')is useful to determine the annihilation parameter A since: A(B'^ — > 
D^K^) = A ) and the average for r£,o/f+ from UTp : 

rrfiK* = 0.10 + 0.04 @ 68% probability 

and we decide to accept, in case of ambiguity, both the solutions (the one with the sign + and 
the one with the sign -), we get : 

roojfo = 0.27 + 0.18 ( @ 68% probability) 

Where, of the error on roo/^o, a contribution of 0.08 is due to the error on vjyiK* and 0.16 to the 
other uncertainties. For the ratio '-^^ : 



= 0.43 + 0.19 (@68%probability) 



If we extrapolate the errors on the measurements to a statistics of 500 fb ^ we would have 
ambiguity in 7.2% of cases. In this situation we would get: 

rjyogo = 0.24 + 0.15 ( @ 68% probability) 

Where, of the error on roogo, a contribution of 0.07 is due to the error on roOK* and 0.13 to the 
other uncertainties. For the ratio : 

= 0.44 ± 0.18 ( @ 68%probability) 

In principle, if one could determine the ratio r^jo^fo, this would be a useful input in the analysis of 
sin(2yS+y). We made a study to see the sensitivity to sin(2yS+7). This study was made assuming 
the information we get on rjyogo in the case of no ambiguity (~ 93% of cases for 500 fb-^),for 
which we get an error on roojfo of 0. 12. Assuming: 

rjyOKO = 0.30 + 0.15 ( @ 68% probability) 

and errors on the observables S and C [55] 

6S = 0.6 6C = 0.5 



we get an error of 60% on the determination of (2/3 + y) through this method. The method 
doesn't seem much sensitive to (2/3 + y): the same exercise made assuming a statistic of 1 ab~^ 
and a value of rjyogo = 0.40 returns an error on (2jS + y) of 42%. 
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Figure 10: y determination with incomplete input {i.e. cases when the number of observables is 
less than the number of unknown parameters). The upper horizontal line corresponds to low- 
luninosity i.e. around current B-factories whereas the lower horizontal curve is relevant for a 
SBF. Blue uses all CPES modes of D^, red is with only K^n~ and purple uses combination of 
the two. Green curve again uses on D^, K*n~ but now includes K*~ and D*^; see text for 

details. 

4 Measurements of 7/^3 

4.1 Combined Strategies for extraction of a clean 7 

Contribution from D. Atwood 

A promising feature of clean and precise extraction of y from B KD type of modes is 
that a multitude of strategies exist which are very effective when used in combination [5 6]. In 
here we focus on the use of direct CP from B "K'"D^ /D^ though it is also possible to use 
time dependent CP violation via 5" "K'^"D^. 

Fig 10 illustrates several of the important features. Here xlmn versus y is compared for 
various methods and input data sets. Let us focus first specifically to B' K'D^, with 
D° decays to common modes such as (CPNES) K^n (ADS[57]) or (CPES) K^^n^ (GLW[58]). 
Then for each such common mode (say K'^n~) there are basically three unknowns: 6st, rg and y 
where 6st is the strong phase and rg = Br{B~ K~D^)/Br{B~ K~D^) and two observables: 
the rate for B" and for B^ . Therefore as such we cannot hope to solve for y. However, as 
we add another common mode of say K*p~ , then one is adding one new unknown (a 

strong phase) but 2 more observables. So now there are 4 observables and 4 unknowns and the 
system, is in principle, soluble though in practice presence of discrete ambiguities complicates 
the solution. In general for N such common modes, there will be N + 2 unknowns and 2N 
observables and as more modes are added y can be solved for very effectively. 

Fig. 10 and Fig. 1 1 show situation with regard to under-determined and over-determined 
cases respectively. The upper horizontal line corresponds roughly to the low luminosity i.e. 
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Figure 1 1 : y extraction with over-determined cases. Purple curve shows the effect of combining 
GLW (all CPES modes) with one ADS (K'^n~) mode; black curve differs from purple only in that 
it also includes from Z)*°; blue curves show the effect of properly including the correlated 
strong phase between D^* + n and Z)"* ^ + y. Orange curve includes lot more input 

including Dalitz and multibody modes, see text for details (See also Fig 10). 

comparable to the current 5-factories[59] whereas the lower horizontal curve is relevant for a 
super 5-factory. In Fig. 10 in blue is shown the case when only the input from (GLW) CPES 
modes of D° is used; note all the CPES modes are included here. You see that the resolution on 
7 then is very poor. In particular, this method is rather ineffective in giving a lower bound; its 
upper bound is better. 

In contrast, a single ADS mode (A'^/r ) is very effective in so far as lower bound is con- 
cerned, but it does not yield an effective upper bound (red). 

Note that in these two cases one has only two observables and 3 unknowns. In purple is 
shown the situation when these two methods are combined Then at least at high luminosity 
there is significant improvement in attaining a tight upper bound; lower bound obtained by 
ADS alone seems largely unaffected. 

Shown in green is another under determined case consisting of the use of a single ADS 
mode, though it includes K* as well D*"; this again dramatically improves the lower bound. 
From an examination of these curves it is easy to see that combining information from different 
methods and modes improves the determination significantly [56]. 

Next we briefly discuss some over determined cases (Fig. 11). In purple all the CPES modes 
of are combined with just one doubly cabibbo-suppressed (CPNES) mode. Here there are 
4 observables for the 4 unknowns and one gets a reasonable solution at least especially for the 
high luminosity case. 

The black curve is different from the purple one in only one respect; the black one also 
includes the D*'* from B~ K D*^* where subsequently the D°* gives rise to a D". Comparison 
of the black one with the purple shows the remarkable improvement by including the Z)"*. In 
this case the number of observables (8) exceeds the number of unknowns (6). 

Actually, the D^* can decay to via two modes: D"* + n or + y. Bondar and 
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Gershon [60]have made a very nice observation that the strong phase for the y emission is 
opposite that of the n emission. Inclusion (blue curves) of both types of emission increases the 
number of observables to 12 with no increase in number of unknowns. So this improves the 
resolving power for y even more. 

The orange curves show the outcome when a lot more input is included; not only K~ , K~*, 
D^, D^* but also Dalitz and multibody decays of are included. But the gains now are very 
modest; those once the number of observables exceeds the number of unknowns by a few (say 
0(3)) further increase in input only has a minimal impact. 

4.2 Measurements of 7 at BaBar 

The measurement of the angle 7 of the unitarity triangle through B~ /)(*)/^(*)- requires the 
combination of as many B and D modes as possible to reduce the statistical uncertainty. Though 
the most statistically powerful approach analyzed so far is the one exploiting the Dalitz plot 
analysis of K^^n'Tj:* from B' D^*^^K-, the GLW and ADS methods also provide useful 
information. 

4.2.1 BaBar GWL and ADS 

Contribution from M. Rama 

BaBar has studied the decays B^ £)(*)(' £>*() _^ £)0^o _^ ^0^- 

Z)° meson is reconstructed in the CKM-favored modes K~n^, K~n^n^n' and K'n'^n^; the CP+ 
eigenstates K~K'^ and n'n'^; and the CP- eigenstates K^n^, K^(p, K^oj. Figure 12 shows the dis- 
tributions of AE of the B —> D^K events (CP+ and CP- modes). The signal yields are extracted 
through a maximum likelihood fit that uses as input AE and the Cherenkov angle of the bache- 
lor track K. Similar techniques are used to select B D*^^K and B D^K* [61, 62, 63, 64]. 
The measurements of Acp± and R± allow to constrain the unknowns r/,, 6b and the CKM an- 
gle 7. The results of the measurements are reported in Table 6, where Rfiav = BF(B~ —> 
D^*^^K~)IBF{B~ D^*^^7T~). Particular care is needed when evaluating the systematic un- 
certainties associated to peculiar sources of background. The B~ D^p^K~ decays are af- 
fected by the charmless 3-body background B~ K~h~h* (h = n, K), characterized by the 
same mes and AE distribution as the signal. The B~ K~n~n* and B~ K K back- 
grounds (4 + 4 and 29 + 7 events, respectively) are subtracted from the B~ D^\-^ n~n*'\K~ 
and B~ > K~K'^]K~ signals (18 + 7 and 75 ± 13 events, respectively). Similarly, the 

B- D%-^ K'^^aoW*- is evaluated and subtracted from the B' D^[-^ K^^cpW signal. In 
this case, due to the different spin-parity of ao with respect to (p, the background has an opposite 
CP asymmetry with respect to the signal. 

(-)„ 

On a datasample of 227 million T(45 ) BB decays we have searched for B D K 

(-) 

followed by K^n , as well as the charge conjugate sequences. In these processes, the 

favored B decay followed by the doubly CKM- suppressed D decay interferes with the sup- 
pressed B decay followed by the CKM-favored D decay. The yields of the signal mode and 
the normalization mode (B~ D^K~ with K~n*) are extracted through a fit on m^s , 

after requiring that the events have a AE value consistent with zero within the resolution. We 
find a total of 4.7;|;3 2 signal events and 356 + 26 normalization events, that are used to evaluate 
the charge-integrated ratio Rsn = [5(5" ^ [K^7t-]dK-) + B(B^ [K-7t^'\dK^)]/(2B(B^ ^ 
[K'^tt^IdK^)) = rl + rl + IrBrj) cos (^b cos 7. The resulting limit is I^Kn < 0.030 at 90% CL. 
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Figure 12: Distributions of IS.E for events enhanced in 5 ^ D^K signal. Left: CP+; right: 
CP-. Solid curves represent projections of the maximum likelihood fit; dashed-dotted, dotted 
and dashed curves represent the B D'^K, B D'W and background contributions. 



Using rjj = 0.060 + 0.003 [66], and allowing the variation of y and 6b on the full range 0° - 180°, 
we set the limit < 0.23 at 90% CL [65]. 

Table 6: Measured CP asymmetries Acp± and ratios R+. The first error is statistical, the second 
is systematic. 



B- D^K- B- D^K*- B D*^K- 



A^B5(xlO^) 214 227 123 

Acp+ 0.40 + 0.15 + 0.08 -0.09 + 0.20 + 0.06 -0.02 + 0.24 + 0.05 

Acp- 0.21 +0.17 + 0.07 -0.33 + 0.34 + 0.10' 1.09 + 0.26!°^° 

i?+ 0.87 + 0.14 + 0.06 1.77 + 0.37 + 0.12 *** 

R 0.80 + 0.14 + 0.08 0.76 + 0.29 + 0.06 1 *** 



5 Measurements of 7/^3 with D" Dalitz analysis 

The most precise determination of the standard model CKM phase in the long run is provided 
by the methods based on the interference between b cus and b ucs [58, 46]. In the case of 
charged B decays this means that the interference is between B' DK~ followed by D ^ / 
decay and B DK' followed by D ^ /, where / is any common final state of D and D 
[58, 46, 60]. Here we will restrict ourselves to the case, where / is a multibody final state. For 
concreteness we focus on the following cascade decay [17] 

B- DR- (Ks7r-n'-)DK-. (31) 

For the B decay amplitudes we define 

A(B-^D'K-) = Ab, (32) 

A(B dV") = ABrBe'^^'-^\ (33) 

'Additional biases 6Acp- = 0.15 ± 0.10 ■ (Acp- - Acp+) and SR- ^Zq^u are quoted by the authors, reflecting 
possible interference eiJects between the and a> resonances and the background. 
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Here 5^ is the difference of strong phases and is taken to be positive. The same definitions 
apply to the amplitudes for the CP conjugate cascade DK^ {Ks n*n~)DK^ , except that 
the weak phase flips the sign 7 ^ -y in (33). For the D meson decay we further define 

AD{sn,sn)=An,ne''''-'' = A(D' ^ Ksipi)n~ip2)n^ip3)) (34) 

= A(D ^ Ks{p,)n^(p2)n-(p3)), 

where Sij = {pi + pjf, and pi,p2,P3 are the momenta of the Ks,n' respectively. Note that 
in the last equation CP symmetry was used. 

Obviously, if A^isn, su) were known, one could extract y from B- DK- decay widths. 
One can measure Ao{si2, ^'13) from tagged D Ksn^n by modeling it with a sum of Breit- 
Wigner forms and fitting the parameters from the corresponding Dalitz plot. This is the ap- 
proach used at present by both BaBar [73] and Belle [74, 75], where the modeling error on 
extracted value of y is estimated to be around 10°. In the future this modeling error can be 
avoided by performing a model independent analysis [17, 72]. 

Once the function A£)(m^, m^) is fixed using a model for the Ksn^n' decay, the Dalitz 
distributions can be fitted simultaneously using the expressions for the two amplitudes and rg, 
6 and y can be obtained. 

To illustrate the region of the Dalitz plot most sensitive to y measurement, we show in 
Fig. 13 the distribution of simulated B~ D^K~ events based on our Dalitz model, where each 
event is given a weight of Here L is the likelihood function described in the following 

section. The regions of interference between doubly Cabibbo suppressed and Cabibbo allowed 
decays and CP eigenstate decays are clearly the most sensitive ones: they exhibit the highest 
weights. 




m^2{GeV^/c'') 



— _ — 

Figure 13: Z) Kgn'^n Dalitz distribution of simulated (B^ D K^) events. Each event is 
given a weight The weight scale is indicated on the right of the plot. The black points 

represent the same events with weight equal to unity. 
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Figure 14: Dalitz distributions from Belle, separately for and B decays, for (left) 5* 
DK-, (middle) B- D*K-, (right) B- DK*-. In each pair, the left plot contains the 5+ 
candidates, while the right shows those from B~. 



5.1 Belle Measurement of 03 From B- -> D^*^K^*^-, D -> Ksn^n' Dalitz 
Analysis 

Contribution from K. Abe 

Belle have measured 03 using the Dalitz analysis method in B- £)(*)^(*)± using a data 
sample of 275 million BB pairs. The subdecays D* DtiP, K*- Ksn-, D Ksn^n' are 
used. The numbers of candidate events are 209 + 16 (with a background fraction of 25 + 2%) 
for 5* ^ DK*, 58 + 8 (13 + 2%) for 5* ^ Z)*is:*, and 36 + 7 (27 + 5%) for 5* ^ D/sT**. The 
Dalitz distributions of candidate events are shown in Fig. 14. 

We use continuum D*- [^^/r^/r Jd/t* decays to obtain a flavour tagged D" — > Ksn^n~ 
sample and express the Ksn*n~ amplitude as a sum of 18 resonant and one non-resonant 
amplitudes. 

N 

f(ml, ml) = aje'"'Aj(ml,ml) + (35) 
i=i 

The expected Z)° Ksn*n~ Dalitz distribution, with which we compare the data, is then given 
by 

\f{ml + n^, ml + ^'■fQxp{-—^-^)dix^ + B(m^, m!) (36) 

with efficiency e, tttt mass resolution cr^, and background B. We obtain the parameters of the 
Dalitz model by minimizing 



-21ogL = -2 



n p 

y ^ log /'(m+ i) - log I p(m^, ml)dm^dml 



(37) 



with free parameters and aj for each resonance (except for Ksp for which aj = 1, or^ = 
are reference values), b, /3, and the masses and widths of cri and cri. The fit results are listed in 
Table 7. 

The fitting procedure for B- £)(*)^*± is similar to that for the D" Kgn^n' fit ex- 
cept f{m\, m^) is replaced with /(m^, ml) + re^'^^'f'^^^^ f{ml, m\). The efficiency, resolution, and 
background distributions are also replaced with those relevant for the B decay. 

We use 5* D^*^n*' as control samples to test the analysis procedures. Separate fits to the 
5+and5- data give r+ = 0.039+0.021,0+ = 240°+28°, and r_ = 0.047+0.018, 6'_ = 193°+24° 
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Intermediate state 


Amplitude 


rhase ( ) 


tit traction 


KsCTx 


1 cn 1 A 1 A 

1.57 ± 0.10 


O 1 /I 1 A 

214 ± 4 


9.0% 


1/- 


l.U (nxed) 


(nxed) 


21.6% 


J/' ^ , 


U.U31U ± U.UUlU 


111/1 1 in 

113.4 ± i.y 


A /I (97 

{jAvc 


A.5/o(y«o) 


u.jy4 ± U.UUO 


OAT 1 '3 


A CiCf 




U.zJ ± U.UJ 


O 1 A 1 1 Q 

ZlU ± 1 J 


A ^^I)/ 
U.D% 


i^5/2(1270) 


1 O 1 A A /I 

1.32 ± 0.04 


O /I O 1 O 

348 ± 2 


1 C/77 

1.5% 


KsMUlO) 


1.25 + 0.10 


69 + 8 


1.1% 


T/- 0/1/1 cfw 

i^sp"( 1450) 


A OA 1 A C\'~l 

0.89 + 0.07 


1 1 £ 

1 + 6 


A /I /V 

0.4% 


A (892) TT 


1 1 1 A A1 A 

1.621 + 0.010 


111 T 1 AC 

131.7 + 0.5 


61.2% 


a: (092) TV 


A 1 C /I 1 A AAC 

0.154 ± 0.005 


'J 1 ^7 ^7 1 1 ^ 

31 /. / + 1.6 


0.55% 


A (1410) ;r 


0.22 + 0.04 


120 + 14 


0.05% 


t^* /" 1 /I 1 A\ — — + 

A (1410) ;r 


A O C 1 A A/1 

0.35 + 0.04 


253 ± 6 


0.14% 


Aq(1430) ;7r 


O 1 C 1 A A/1 

2.15 + 0.04 


n AO '1 1 1 1 

348.7 + 1.1 


T /I /7/ 

7.4% 


Aq(1430) TV 


A CO 1 A A/I 

0.52 + 0.04 


OA 1 A 

89 + 4 


A A'^frJ 

0.43% 


A2(1430) ;r 


111 1 A AO 

1.11 + 0.03 


O O A C 1 1 O 

320.5 + 1.8 


2.2% 


v*('\ /l'?^^~^r+ 

Aj^l'+JUJ TT 


U.ZJ ± u.uz 


OA'? -1- 1 
Z.OJ + / 




i^*(1680)+;r- 


2.34 + 0.26 


110 + 5 


0.36% 


if*(1680)-;r+ 


1.3 + 0.2 


87+11 


0.11% 


non-resonant 


3.8 ± 0.3 


157 ±4 


9.7% 



Table 7: Results of the fit to obtain the parameters of the Z)° Ksit^n decay model, from 
Belle. 



for the T>n- samples, and r+ = 0.015 + 0.042, 6i+ = 169° + 186, and r_ = 0.086 + 0.049, 
Q- = 280° + 30° for the D*n- samples. These results are not inconsistent with the expectation 
of r ~ \VubV*^/VcbV*J ~ 0.01 - 0.02 although ~ 2cr deviation is possible, which we include in 
the systematic error. 

Fig. 15 shows the results of the fits in terms of Re(re'^)-Im(re'^) separate B'^ and B~ fits for 
DA'*, D*K-, and DK*- samples, respectively. 

We use the frequentist technique for determining the confidence region. Here we determine 
PDFs using a toy Monte Carlo and calculate the confidence level and estimate the true (r, (pj, 6) 
and their errors. This procedure removes possible bias which arises from an assumption of r 
being positive-definite. For the 5 and ^3, the diff'erence from unbinned fits are small. However, 
r is changed by a small amount as expected. The results are shown in Table. 8. 

Systematic errors arise from, i) background fractions and their shapes in the Dalitz plane, 
ii) efficiency shapes, iii) momentum resolution shapes, and iv) possible bias in the analysis 
method. In addition, another error is assigned due to the model-dependence of the variation 
of the complex phase in the —> Kstt^tt' amplitude. Furthermore, for the DK** mode, 
we include an additional systematic error due to the possible presence of non-resonant DKstt 
component which can give interference with different values of r and 6. 

When the DK- and D*K- modes are combined, we obtain ^3 = 68°!j5l + 13° + 11°. In 
the near future, the statistical error can be improved by combining all modes, DAT*, D*(Dn)K'^, 
DK*"*, D*{Dy)K'*, and by exploring other D decay modes such as D^ n*n~n^ and D*^ 
K^K'Ks. As the data size increases, the experimental systematic error, which is presently 
dominated by the D^*^n control sample size, will also improve. We hope to improving the D° 
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Figure 15: Results from Belle of separate 5+ and B fits in terms of Re(re'^^)-Im(re'^) for (left) 
5± ^ DK^, (middle) 5± ^ D*K^, (right) 5± ^ DK*^. CP violation is expected to appear as a 
relative rotation around the origin of the B'^ and B~ contours. 





r 


03 (°) 


5(°) 




0.21 + 0.08 + 0.03 + 0.04 
0.12!°j^ + 0.02 + 0.04 
0.25:" }^ + 0.09 + 0.04 ± 0.08 


64+19+13 + 11 
75 + 57 + 11 + 11 
112 + 35 + 9 + 11 + 8 


157 + 19 + 11+21 
321 +57 + 11 + 21 
353 + 35 + 8 + 21 +49 



Table 8: Results of the Belle 5± d^*^K^*^^, D ^ Ksn^n' measurements. 



decay model by studying the CP-tagged D Ksn^n' data from (^(3770) D°D^ at CLEO-c 
or BES. These data will provide |/(m+,m^) + f{mi,m\)f' and therefore a direct measurement 
of the phases, thus dramatically reducing the model uncertainty. 

5.2 Extracting 7 from D Dalitz analysis at BaBar 

Contribution from M.H. Schune 

5.2.1 The Dalitz Model 

The function Ad has been obtained using the flavour tagged D meson sample from the contin- 
uum decays D*- Dn^. The Dalitz (ml, ml) distribution (Fig. 16) is fitted in the context of 
the isobar formalism described in [76]. In this formalism the amplitude / can be written as a 
sum of two-body decay matrix elements and a non-resonant term according to the expression: 

= ^nre"^"' + Z^are"^' A,{Ksn' n^\r). Each term of the sum is parameterized with an amplitude 
and a phase. We fit the Dalitz distribution with a model consisting of 13 resonances leading 
to 16 two-body decay amplitudes and phases (Table 9). Of the 13 resonances eight involve a 
Ks plus a nn resonance and the remaining five are made of a {Ksn~) resonance plus a n'^. We 
also include the corresponding doubly Cabibbo-suppressed amplitudes for most of the {Ksn~) 
n'^ decays. All the resonances considered in this model are well established except for the two 
scalar nn resonances, cti and 0-2, whose masses and widths are obtained from our sample. 

An unbinned maximum likelihood fit is performed to measure the amplitudes a„^, and the 
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Figure 16: (a) Ksn n'^ Dalitz distribution from D** D^n'^ events, and projections on 

(b) m^, (c) ml, and (d) m^„+„- \ from BABAR. The curves are the fit projections. 



phases (pnr, (f>r- components. The results of the fit are shown in Fig. 16. Amplitudes, phases and 
fit fractions as obtained by the likelihood fit are reported in Table 9. 



Resonance 


Amplitude 


Phase (deg) 


Fit fraction 


if* (892)- 


1.781 


+ 0.018 


131.0 + 0.8 


0.586 


k;(1430)- 


2.447 


+ 0.076 


-8.3 + 


2.5 


0.083 




1.054 


+ 0.056 


- 54.3 + 


2.6 


0.027 


is:*(i4i0)- 


0.515 


+ 0.087 


154 + 


20 


0.004 


K*il680y 


0.89 


+ 0.30 


- 139 + 


14 


0.003 


K*(S92y 


0.1796 


+ 0.0079 


-44.1 + 


2.5 


0.006 


is:,*)(1430)+ 


0.368 


+ 0.071 


-342 + 


8.5 


0.002 




0.075 


+ 0.038 


- 104 + 


23 


0.000 


P(770) 


1 (fixed) 


(fixed) 


0.224 


co(lS2) 


0.0391 


+ 0.0016 


115.3 + 


2.5 


0.006 


/o(980) 


0.4817 


+ 0.012 


-141.8 + 


2.2 


0.061 


/o(1370) 


2.25 


+ 0.30 


113.2 + 


3.7 


0.032 


/2(1270) 


0.922 


+ 0.041 


-21.3 + 


3.1 


0.030 


p(1450) 


0.516 


+ 0.092 


38 + 


13 


0.002 


cr 


1.358 


+ 0.050 


-177.9 + 


2.7 


0.093 


a' 


0.340 


+ 0.026 


153.0 + 


3.8 


0.013 


Non Resonant 


3.53 


+ 0.44 


127.6 + 


6.4 


0.073 



Table 9: Amplitudes, phases and fit fractions as obtained by the likelihood fit on the tagged D 
sample, from BABAR. 

We estimate the goodness of fit through ax^ test and obtain = 3824 for 3054-32 degrees 
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Figure 17: Contours at 68.3% (dark) and 95% (light) confidence level (statistical only) in the 
{x^*)^y(*)') planes for (a) DK~ and (b) D*K~, separately for B~ (thick and solid) and (thin and 
dotted). Two-dimensional projections in the r^*' - y planes of the five-dimensional one- (dark) 
and two- (light) standard deviation regions, for (c) DK~ and (d) D*K' . 



of freedom. 

5.2.2 Results of the CP fit 

BaBar reconstruct the decays B~ DK~ and B' D*K~ with D* Dn^ , Dy. B~ candi- 
dates are formed by combining a mass-constrained D^*^ candidate with a track identified as a 
kaon [33]. 

We simultaneously fit the B~ D^*^K~ samples using an unbinned extended maximum- 
likelihood fit to extract the CP-violating parameters along with the signal and background 
yields. 

In the sample of 21 1 million BB events, the following signal yields are obtained 

NiB~ D'^K~) = 261 + 19, (38) 
N(B- D*"(DV')i^") = 83 + 11, 
N(B~ D*^{D\)K-) = 40 + 8. 

The results for the CP-violating parameters and y+\ are defined as the real and imagi- 
nary parts of the complex amplitude ratios r^*^e'^''s'*^\ respectively. 

The and variables are more suitable fit parameters than r^*\ and y because they 
are better behaved near the origin, especially in low- statistics samples. Figures 17(a,b) show the 
68.3% and 95% confidence-level contours (statistical only) in the {x^*\y^*^) planes for DK' and 
D*K~, and separately for B and B^. The separation between the B~ and B^ regions in these 
planes is an indication of direct CP violation. 

The largest contribution to the systematic uncertainties in the CP parameters comes from 
the choice of the model used to describe the D° Ksn'n'^ decay amplitudes. To evaluate this 
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uncertainty we use the nominal Dalitz model (Table 9) to generate large samples of experiments 
and we compare experiment by experiment the values of and obtained from fits using the 
nominal model and a set of alternative models. Models where one or both of the cr resonances 
are removed lead to a significant increase in the of the fit. We use the average variations of 
x^*^ and y^!^ corresponding to this second set of alternative models as the systematic uncertainty 
due to imperfect knowledge of A/). 

A frequentist construction of the confidence regions of p = {rB,r*g, 63,8*^,7) based on the 
constraints on x^*^ and _y±^ has been adopted (see below and Ref. [66]) and allows to extract 7, 
vb, 4, 6b and (5^: 

7 = (70 ± 31 !J|)° [12°, 137°], re = 0.118 + 0.079 + 0.034 [0,0.279], r% = 0.169 + 

0.096 toms -aoi [0, 0.354], 6b = (l04 + 45 ^b = (296 + 41 1}^ + 15)°. The first error 

is statistical, the second is the experimental systematic uncertainty and the third reflects the 
Dalitz model uncertainty. The values inside square brackets indicate the two-standard-deviation 
intervals. 



5.3 7 extraction in B- -» D^*^K- Dalitz: a frequentist statistical treat- 
ment 



Contribution from N. Neri 



5.3.1 Cartesian coordinates as fit parameters 



The likelihood function used for the measurement of the angle 7, reconstructing D° three-body 
decays (such as D° Ksn^n~) from B~ —> D^*^K~ decays through a Dalitz analysis [17], is 
affected by biased values and non-Gaussian errors if we fit directly for the parameters rg, 7 and 
6b- We find that cartesian coordinates: 

x± = rB cos(7 + 6b) 
= rB sin(7 + 6b) 

are more suitable set of parameters because they are better behaved near the origin expecially 
in low statistics sample. In addition the cartesian coordinates are largely uncorrected. If we 
represent in the ix±,y±) plane the results for B- decays, the distance of the two points d is 

d = [(x_ - + (y_- y^fyn = 2rB\ sinj] . 
A non null distance means evidence of direct CP violation. 

5.3.2 Frequentist procedure 

In the case of a single signal sample^ DK"^, we define a 4-dimensional PDF, of the fit 

parameter as a function of the true parameter, as: 

^5^-^(z+,z_|p*) = G2{z+;r'BCOs(6' + y),r'B^HS' + y),o-^,,o-y^,p+) x 

G2 (z_; cos((5' - Y), 4 sm(6' - y'), cr;,_, cry_,p_) 



'a similar procedure has been used for combined measurements of two and tliree signal B"^ — > 
samples. 
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where 



1 


^^-^tx•r 


2(l-p2) 


\ 4 



and z+ = >'±) P = (^b, 7, S). The vectors z+ and defined equivalently to z+ and p re- 
spectively, are the corresponding parameters in the truth parameter space. The Gaussian widths 
{cTx^, CTyJ and the correlation distributions (p±) can obtained either from Toy MC experiments 
or from the fit to the data sample itself. 

The confidence level 1 - or for each set of true parameters is calculated as 



where the integration domain D (the confidence region) is given by the condition 
-^^(z,,z_|pO > ^^^(zf «,z1-1pO , 

i.e. it includes all points in the fit parameter space closer to the truth point than the data point 
^dato fffQ are interested in building a 3-dimensional region of joint probability 1 - ao, then we 
select only those points for which aip') < a^. The 2-dimensional and 1 -dimensional contours 
are then built by projecting the 3-dimensional joint probability regions. The correct coverage of 
the method has been verified. 

Finally, the significance of CP violation can be determined by finding the confidence level 
1 - acp for the most probable CP conserving point p^, i.e. the p' point with = or y = 0. 

5.3.3 Perspectives 

A preliminary study has been performed for the —> K^K^K' decay mode. The B~ candidate 
is reconstructed using similar selection of the KsTt'tt'^ mode Using this event yield and a 

study of the KsK'K^ Dalitz plot, toy studies indicate a marginal gain of few percent on 

the 7 uncertainty. However, it is worthwhile noting that the systematic uncertainty due to the 
Dalitz model will be totally independent of the one of the —> Ksn~n^ analysis. 

A first analysis of the B~ D^K~ mode using the —> n^n~7fi decay mode has been 
performed. This channel is aff'ected by a larger background that could dilute the sensitivity to 
y. Nevertheless, adding this channel is certainly worthwhile in a strategy of a combined fit with 
other channels. 

Using the present cartesian central values and a Bayesian technique we have computed the 
expected uncertainty for y for various integrated luminosities. The systematic uncertainty has 
been (conservatively) assumed to stay constant. It has been shown that, up to 500 fb~\ the 
result is yet not dominated by the Dalitz model systematic uncertainty. However, in a view 
of a larger statistics accumulated by 5-factory it is probably worthwhile performing a model 
independent analysis [17]. 

5.4 Model independent approach 

Contribution from J. Zupan 
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The first step in the model independent approach is to partition the Dalitz plot into 2k bins 
placed symmetrically about the 12 13 line. The k bins lying below the symmetry axis are 
denoted by index i, while the remaining bins are indexed with i. The i-th bin is obtained by 
mirroring the i-th bin over the axis of symmetry. For bins below the symmetry axis we define 

Ci = ^i/p Ai2,i3 ^13,12 COS((5i2,l3 - <^13,12), (39) 

Si = JdpAi2,i3M3,nsm(6i2,u-Si3,n), (40) 
Ti ^ JdpAl^,,, (41) 

where the integrals are done over the phase space of the i-th bin. The variables c, and 5, contain 
difi'erences of strong phases and are therefore unknowns in the analysis. The variables r„ on 
the other hand, can be measured from the flavor tagged D decays, and are assumed to be known 
inputs into the analysis. The variables c„ 5/ of the i-th bin are related to the variables of the i-th 
bin by 

cj = Ci, S] = -Si, (42) 

while there is no relation between T, and Tj. 

Together with the information available from the B"^ decay, we arrive at a set of 4k equations 

tj = Jdr(B- ^ (Ksn-n*)DK-) = 

Ti + rlTj + 2rB[cos(5B - y)ci + sm{6B - y)si\, (43) 

where we display only the first k equations, while the other expressions for f7, f^^, can be 
obtained from the above by 12 13 and/or y -y replacements [17]. We have 2k + 3 
unknowns: c,-, Si, rg, 6b, J for Ak observables, so that the set of equations is solvable for k>2, 
i.e., if Dalitz plot is divided into at least four bins. The whole approach has been extended also 
to neutral B decays [77]. 

An important input can be provided by charm factories [17, 72]. Namely, the parameters 
Ci and s^^ describing D decay can be measured at charm factories working at il/(3nO). This 
greatly reduces the number of unknowns needed to be fit from B decays (which is the limiting 
statistical factor). Another observation is that c,, are bounded from above and below, I*;!, \ci\ < 
J.dp Ai2,i3Ai342 < V^/ Ti- This can prove useful in the actual implementation of the method. 
An important question for the implementation of the method is also, how small need the bins be, 
not to average out sensitivity to y? Answer to this question is rather complicated and difficult 
to answer without a Monte Carlo study. 

Another interesting observation is, that to leading order in x = Amo/F/), y = AF£)/(2F£)) 
CP conserving D-D mixing does not change the formalism [78]. Also, the uncertainty due to 
CP violation in D sector is ~ 10""^ suppressed in SM and therefore completely negligible. 
However, even these eff'ects can be included in the analysis [78]. 

5.5 Charm phenomenology and extraction of the CKM angle 7 

Contribution from A. A. Petrov 
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The Standard Model is a very constrained system, which implements a remarkably simple 
and economic description of all CP-violating processes in the flavor sector by a single CP- 
violating parameter, the phase of the CKM matrix. This fact relates all CP-violating observables 
in bottom, charm and strange systems and provides an excellent opportunity for searches of 
physics beyond the Standard Model. One of the major goals of the contemporary experimental 
B physics program is an accurate determination of the CKM parameters. As we shall see below, 
inputs from charm decays are important ingredient in this program, both in the extraction of the 
angles and sides of the CKM unitarity triangle. 

The cleanest methods of the determination of the CKM phase 7 = arg [-^Md^*^/^cci^*^] 
involve the interference of the b —> cus and b —> ucs quark-level transitions [79]. A way to 
arrange for an interference of those seemingly difi'erent processes was first pointed out in [58, 
46]. It involves interference of two hadronic decays 5"^ D^K^ fK* and D^K^ 
fK*, with / being any common final state for and decays. Since then, many different 
methods have been proposed, mainly differing by the fK final state (i.e. with / being a CP- 
eigenstate or not a CP-eigenstate) and paths of reaching it: a combination of the Cabibbo- 
favored (CF) and doubly-Cabibbo suppressed decays (DCSD) D°(D°) ^ K^n' [57], singly- 
Cabibbo suppressed decays (SCSD) D*'(Z)") KK* [67], Cabibbo-favored decays employing 
large K*^ - mixing transitions [17]. For similar methods involving interference of the initial 
state, see [80]. 

All these methods produce expressions that depend on y and several hadronic parameters, 
such as ratios of hadronic amplitudes and strong phases. These parameters cannot be reliably 
computed at the moment, so their values must be fixed from experimental data. In principle, 
all of the hadronic parameters in these methods can be obtained from the measurements of 
5-decays only. However, accuracy of these methods can be significantly improved if some 
measurements of charm-related parameters are performed separately. A good example is pro- 
vided by the original GW method [58, 46], which does not take into account the possibility 
of relatively large - mixing [81]. This method relies on the simple triangle amplitude 
relation 

V2A(5+ ^ D^K^) = A(B^ + A(5+ ^ Wk^), (44) 

which follows from the relation which neglects CP-violation in charm, 

^f2\D^) = + \D^). (45) 

An amplitude A(B* D+K*) is measured with D decaying to a particular CP-eigenstate. 
Neglecting - mixing, angle 7 can then be extracted, up to a discrete ambiguity, from the 
measurements of 5* —> fcpK"^ and 5* —> D^, D^K"^. In particular, 

r[5* ^ fcpm oc I +rl + Imc^, (46) 

where c± = cos(7 + 6b) and and are defined from 

A{B* D^K^)/A(B^ Wk*) = ree'^^^'^'K (47) 

Then, 

1 r n — m ^^^^ 



• 2 1 
sm 7 = - 

' 2 



l-c,c_± ^(l-ci)(l-c!) 



It is easy to see that - mixing, if not properly accounted for, can affect the results of 
this analysis. Indeed, taking - mixing into account results in the modification of the 
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definitions of c±. 



cos(7 + 6b) + sin 29d - I'^^ffo cos (y + 26 d + 6b) + cos 26 d \ , 
Ltb /tt^ l j 



2 

2rD 



where 



m2 - mi T2- Fi 

r 2r 



(49) 



(50) 



with mi,2 and Fi 2 being the masses and widths of D-meson mass eigenstates D12 (which reduce 
to Z)± in the CP-invariance limit), 77^ is a CP-parity of fcp, and do is a CP- violating phase of 

- mixing. It is easy to see that y ~ 1% can impact the determination of 7 from these 
modes. Thus, separately constraining - mixing parameters will be helpful. 

The current experimental upper bounds on x and y are on the order of a few times 10^, and 
are expected to improve significantly in the coming years. As was recently shown [82], in the 
Standard Model, x and y are generated only at second order in SU(3)f breaking. 



x,y~ sin^ 6c x [S U(3) breaking]^ 



(51) 



where 6c is the Cabibbo angle. Therefore, predicting the Standard Model values of x and y 
depends crucially on estimating the size of SU(3)/r breaking. Although y is expected to be 
determined by the Standard Model processes, its value can affect the sensitivity to new physics 
of experimental analyses of D mixing [83]. 

Presently, experimental information about the - Z)° mixing parameters x and y comes 
from the time-dependent analyses that can roughly be divided into two categories. First, one 
can look at the time dependence of D —> f decays, where / is the final state that tags the 
flavor of the decayed meson. The most popular is the non-leptonic doubly Cabibbo suppressed 
decay D° K^n~ . Time-dependent studies allow one to separate the DCSD from the mixing 
contribution Z)° ^ ^ K*n~, 



T[D\t) ^ isTV-] = e-^'\AK-nA^ 

X 



Rl 



R + ^fRRmiy' cos - x' sin 0)F? + -j-(y^ + x^XYtf 



(52) 



where R is the ratio of DCS and Cabibbo favored (CF) decay rates. Since x and y are small, the 
best constraint comes from the linear terms in t that are also linear in x and y. A direct extraction 
of X and y from Eq. (52) is not possible due to unknown relative strong phase 6d of DCS and 
CF amplitudes [84], as x' = xcos 6d + ysindo, y' = y cos 6d - xsin^o. As discussed below, 
this phase can be measured independently [81, 85, 86]. Second, mixing can be measured by 
comparing the lifetimes extracted from the analysis of D decays into the CP-even and CP-odd 
final states. This study is also sensitive to a linear function of y via 



-\ = y cos - X sin 



Rl-l 



(53) 



Time-integrated studies of the semileptonic transitions are sensitive to the quadratic form x^ +y^ 
and are not competitive with the analyses discussed above. 

Charm factories (ChF) such as CLEO-c and BES-III introduce new time-independent meth- 
ods that are sensitive to a linear function of y. There, one can use the fact that heavy meson 
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pairs produced in the decays of heavy quarkonium resonances have the useful property that the 
two mesons are in the CP-correlated states [87], 

\DD'')l = l\D%,)D\k2)) + (-lf\D\k2)D\h)}] , (54) 

where L is the relative angular momentum of two D mesons, CP properties of the final states 
produced in the decay of i/^(3770) are anti-correlated, one D state decayed into the final state 
with definite CP properties immediately identifies or tags CP properties of the state "on the 
other side." That is to say, if one state decayed into, say n^Ks with CP = -1, the other state is 
"CP-tagged" as being in the CP = -i-l state. By tagging one of the mesons as a CP eigenstate, 
a lifetime difference may be determined by measuring the leptonic branching ratio of the other 
meson. Its semileptonic width should be independent of the CP quantum number since it is 
flavor specific, yet its branching ratio will be inversely proportional to the total width of that 
meson. Since we know whether this D(k2) state is tagged as a (CP-eigenstate) D± from the 
decay of D(ki) to a final state So- of definite CP-parity cr = +, we can easily determine y in 
terms of the semileptonic branching ratios of This can be expressed simply by introducing 
the ratio 

j^L ^ nilfj^^{H^S^){H^Xl^v)-\ 
" r[ilrL^(H^S^)(H^X)]Br(H^^Xlvy ^ ^ 

where X in H —> X stands for an inclusive set of all final states. A deviation from = I 
implies a lifetime difference. Keeping only the leading (linear) contributions due to mixing, y 
can be extracted from this experimentally obtained quantity, 

3;cos0 = (-l)V:%^. (56) 

Theoretical predictions of x and y within and beyond the Standard Model span several orders 
of magnitude [88]. Roughly, there are two approaches, neither of which give very reliable 
results because is in some sense intermediate between heavy and light. The "inclusive" 
approach is based on the operator product expansion (OPE). In the » A limit, where A is a 
scale characteristic of the strong interactions, AM and AF can be expanded in terms of matrix 
elements of local operators [89]. Such calculations yield x,y < 10"^. The use of the OPE 
relies on local quark-hadron duality, and on A/nic being small enough to allow a truncation 
of the series after the first few terms. The charm mass may not be large enough for these to 
be good approximations, especially for nonleptonic D decays. An observation of y of order 
10"^ could be ascribed to a breakdown of the OPE or of duality. The "exclusive" approach 
sums over intermediate hadronic states, which may be modeled or fit to experimental data [90]. 
Since there are cancellations between states within a given SU(3) multiplet, one needs to know 
the contribution of each state with high precision. However, the D is not light enough that its 
decays are dominated by a few final states. While most studies find x,y < 10 Refs. [90] 
obtain x and y at the 10"^ level by arguing that SU(3);7 violation is of order unity, but the source 
of the large SU(3)/r breaking is not made explicit. It was shown that phase space effects alone 
provide enough SU(3)/r violation to induce y ~ 10"^ [82]. Large effects in y appear for decays 
close to D threshold, where an analytic expansion in SU(3)f violation is no longer possible. 
Thus, theoretical calculations of x and y are quite uncertain, and the values near the current 
experimental bounds cannot be ruled out. 

Another example of a method that will benefit from the separately performed charm decay 
measurements is the ADS method [57]. In its original form, it seeks to obtain y from two sets of 
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measurements, 5* —> D^(—> K^n'^)K'^, i.e. use decays of D into non-CP-eigenstate final states. 



R 



ADS - 



IADS - 



Y{B- ^ Dl^.^K-) + r(5+ ^ Dl^^,-K^) 
To + 2rgrz) cos y cos (Jg + 6d) + r\, (57) 

r(g- ^ Dl„.,-K-) - rjB- ^ D^K-n^K^) 
r(B- ^ + r(5+ ^ Dlj^^^K*) 

= ItbTd cos y cos (5b + 6d) IRads , (58) 

where one parameterizes the ratio of a DCS to CF decays, A(D'' — > K'^7r~)/A(D^ —> K~n*) = 
roe'^'^ . The accuracy of this method will be greatly improved if and So are measured sep- 
arately. While the separate measurement of is already available, determination of 6d will 
became possible at the charm factories at Cornell and Beijing. This allows one to measure 
cos 6d. In order to see this, let us write a triangle relation, 

V2A(D± ^ K'n^) = A(D° K'n*) + A(D° ^ K'n^ (59) 

which follows from the fact that, in the absence of CP-violation in charm, mass eigenstates of 
the neutral D meson coincide with its CP-eigenstates, as in Eq. (45). This implies a relation for 
the branching ratios, 

l±2cosc5.V^ = 2^^1^^^^, (60) 

where we used the fact that yjr^ and neglected CP violation in mixing, which could 

undermine the CP-tagging procedure by splitting the CP-tagged state on one side into a linear 
combination of CP-even and CP-odd states. Its effect, however, is completely negligible here. 
Notice that the phase 6d in Eq. (60) is the same as the one that appears in Eq. (52). Now, if both 
decays of D+ and D_ are measured, cos 6d can be obtained from the asymmetry 

Br{D^ ^ K-n^) - Br(D_ ^ R-jr^) 
cosod = — . (61) 

Both Eqs. (60) and (61) can be used to extract 6d at ChF. Similar measurements are possible for 
other D decays [91]. 

One potential problem in measuring the strong phase this way is the need for high-statistics 

measurement which might not be possible at CLEO-c, provided if phase turns out to be small. 
One indication of that is the fact that A(D" ^ K+n-) = A^A(D^ K^n') = A^A{D K'n^) 
in the flavor SU(3) limit, which implies that (5^ = in this limit. Therefore, calculation of the 
value of 5d is equivalent to computation of an SU(3)-breaking correction. This, however, does 
not imply that it is very small [82], as SU(3) breaking in charm decays can be significant. 

Another possibility would be to use other common modes of D° and D^, such as D ^ 
K*'^7r' ,p'^7T' , etc. The same arguments leading to Eq. (61) can be applied there as well 
and the resulting equation would look identical. One added benefit is that SU(3) symmetry 
arguments alone do not force the strong phase to be zero, so its value could in principle be larger. 
One can speculate that if the chiral symmetry in QCD is realized in the "vector" mode [92], than 
A(D° ^ K*^7r-) = A^A(D" K*-n^) in the limit of "vector" SU(3)lX SU(3)r symmetry, so the 
relevant strong phase is still zero. It is however not clear if this realization of chiral symmetry 
is relevant to charm decays (except, may be in the large A', limit). In any case, vector symmetry 
is badly broken in D decays [92], so the resulting strong phase can still be large. Clearly, only 
experimental measurements would settle these issues. 
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5.6 CLEO-c Impact on 7/ 03 Measurements 



Contribution from D. Asner 

The CLEO-c research program [93] consists of studies of hadronic, semileptonic and lep- 
tonic charm meson decay which include measurements of doubly-Cabibbo processes, studies of 
charm Dalitz plot analyses, and tests for physics beyond the Standard Model including searches 
for charm mixing. Precision determination of 7/^3 depends upon constraints on charm mixing 
amplitudes, measurements of doubly-Cabibbo suppressed amplitudes and relative phases, and 
studies of charm Dalitz plots tagged by flavor or CP content. 

The CLEO-c physics program [93] includes a variety of measurements that will improve 
the determination of 7/^3 from the B-factory experiments. The total number of charm mesons 
accumulated at CLEO-c will be much smaller than the samples already accumulated by the B- 
factories. However, the quantum correlations in the il/{2)llQ)) DD system provide a unique 
laboratory in which to study charm. 

5.6.1 Components of CLEO-c Physics Program Pertinent to 7/03 

Neutral flavor oscillation in the D meson system is highly suppressed within the Standard 
Model. The time evolution of a particle produced as a Z)° or Lf, in the limit of CP conser- 
vation, is governed by four parameters: x = Am/F and y = AT/IT which characterize the 
mixing matrix, 6 the relative strong phase between Cabibbo favored (CF) and doubly-Cabibbo 
suppressed (DCS) amplitudes and Ro the DCS decay rate relative to the CF decay rate. The 
mixing rate Rm is defined as + y^) [94]. Standard Model based predictions for x and y, as 
well as a variety of non- Standard Model expectations, span several orders of magnitude [95]. 
It is reasonable to expect that x » y » 10"^ in the Standard Model. The mass and width dif- 
ferences X and y can be measured in a variety of ways. The most precise limits are obtained by 
exploiting the time-dependence of D decays [94]. Time-dependent analyses are not feasible at 
CLEO-c; however, the quantum-coherent D^D^ state provides time-integrated sensitivity to x, y 
at 0(1%) level and cos 5 ~ 0.1 [93, 85] in 1 fb"^ of (/^(3770) ^ DD. Although CLEO-c does not 
have sufficient sensitivity to observe Standard Model charm mixing the projected sensitivity in 
1 fb~^ at ^^(3770) compares favorably with current experimental results; see Fig. 1 in Ref. [94]. 

5.6.2 Targeted Analyses 

Measurements of each of the four parameters that describe CP conserving mixing can in princi- 
ple be determined individually by a series of "targeted" analyses. The techniques and projected 
sensitivities in 1 fb~^ are described briefly in the subsections that follow. Greater detail can be 
found in Ref. [93]. 

Rm The measurement of Rm can be determined unambiguously by considering the decays 
^/^(3770) ^ (K-n-'XK-n*) and (A(3770) ^ i£^iKX))i£^iKX)). The hadronic final state cannot 
be produced from DCS decay due to a quantum statistics argument - the C-odd initial state can- 
not produce the symmetric final state required by Bose statistics if both the D° and D^ decay 
into the same final state (Kn). A fit utilizing six constraints make this channel efl'ectively back- 
ground free. The double semileptonic channel has fewer constraints due to the two neutrinos. 
The number of "right-sign" (K-n^)(K+n-) and (£-(KX))(e^(KX)) events in the two channels 
combined produced in 1 fb~^ is expected to be ~ 20,000 corresponding to ^2Rm < 1.7% 
@95% Confidence Level (C.L.). 
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COS J Final states where one D decays to Kn and the other to a CP eigenstate can be used to 
probe cos 5. The CP-even eigenstates accessible to CLEO-c include K^K , n^n~ and K'^n^n^ 
and the CP-odd eigenstates include t^n^, K'^i], I^oj, K'^<p, and K'^J^'- The CP content of the 
K^n'^Tr' and n'^n'T^ Dalitz plots will also be utilized. Sensitivity to cos 5 ~l/(2 y/RoN) where 
A^, the total number of CP tagged Kn, will be ~9200 and Rd is already well measured, thus 
leading to an expected precision of +0.09 in cos 6. 

y Tagging one of the D mesons as a CP eigenstate, y can be determined by measuring the 

flavor specific branching ratios of the other meson. The flavor tag width is independent of the 
CP quantum number; however the branching ratio is inversely proportional to the total width. 
Consequently, charm threshold experiments have time-integrated sensitivity to y. Neglecting 
factors related to DCS decays (described in detail in Ref. [96]), we can construct the ratio 

2r[^^(O^Cf)(D^flavor)] 

r[(A^(D^CP)(D^X)5r(D^ flavor)] ^ ^ ^ 

where X represents an inclusive set of all states. A positive y would make the above ratio > 1 
for CP+ tags. 

The decay —> I^n^n' is measured with a Dalitz plot analysis to proceed through inter- 
mediate states that are CP+ eigenstates, such as K^fo, CP- such as K^^p and flavor eigenstates 
such as K*~7:'^ [97]. The presence of mixing through y would introduce an intensity modulation 
across the Dalitz plot as a function of the CP of the contributing amplitudes. The sensitivity at 
CLEO-c to y with Dalitz plot analyses has not yet been fully evaluated. Preliminary estimates 
suggests a limit ofy ~ few% @95% C.L. is attainable with the CLEO-c data. 

Separate measurement of y and Rm as outlined above allow the value x to be determined. 
The upper limit on Rm corresponds to a limit of \x\ < 1.7% @95% C.L.. 

Rp Although Rd is determined by the ratio of {K-n^){K-£^v) to {K' n^){K^ £~v) up to cor- 
rections that are second order in x and y, a more precise measurement is attainable using a D*'^ 
tag to measure D° -> K'^n~ relative to D° -> K~n* [94]. 



5.6.3 Comprehensive Analysis 

The comprehensive analysis simultaneously determines mixing and doubly-Cabibbo suppressed 
parameters by examining various single tag and double tag rates. Due to quantum correlations 
in the C = -1 and C = +1 pairs produced in the reactions e^e' D^^D^\n''^) and 

e^e~ —> D^D^y, respectively, the time-integrated D'^D*' decay rates are sensitive to interference 
between amplitudes for indistinguishable final states. This size of this interference is governed 
by the relevant amplitude ratios and can include contributions from D^-D^ mixing. 
We consider in the comprehensive analysis the following categories of final states: 

• / or /: hadronic states that can be reached from either D*^ or decay but that are not 
CP eigenstates. An example is K~n'^, which results from Cabibbo-favored transi- 
tions or doubly Cabibbo- suppressed transitions. We include in this category Cabibbo- 
suppressed transitions as well as self-conjugate final states of mixed CP, such as K^n'^n'. 

• or €~ : semileptonic or purely leptonic final states, which, in the absence of mixing, tag 
unambiguously the flavor of the parent D. 

• 5+ or 5-: CP-even and CP-odd eigenstates, respectively. 
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Table 10: Estimated uncertainties (statistical and systematic, respectively) for different C con- 
figurations, with rj) and branching fractions to CP eigenstates constrained to the world averages. 



Parameter 


Value 




3 X 10" 




3 X 10" 


pC'=-i _ in . rC=+i _ o y. ia6 


y 





+0.015 


+ 0.008 


+0.007 


+ 0.003 


±0.012 ± 0.005 


x2 





+0.0006 


+ 0.0006 


+0.0003 


+ 0.0003 


±0.0006 ± 0.0006 


cos 6 


1 


+0.15 


+ 0.04 


+0.13 


+ 0.05 


±0.13 ±0.03 


xsinS 









+0.010 


+ 0.003 


±0.024 ± 0.005 



We calculate decay rates for D^'D" pairs to all possible combinations of the above categories of 
final states in Ref. [96], for both C = -1 and C = +1, reproducing the work of Refs. [85, 87]. 
Such combinations, where both D final states are specified, are referred to as double 
tags (or DT). In addition, we calculate rates for single tags (or ST), where either the or 
is identified and the other neutral D decays generically. Any DT event also contains two ST 
decays, and ST rates are obtained by summing the corresponding DT rates. 

Experimental measurements of branching fractions rely on determining yields of ST 
decays assuming knowledge of the luminosity and D^D^ cross section. If, in addition, one also 
measures DT yields, then the luminosity and D^D^ cross section need not be known [98, 99, 
100]. However, in correlated D^D^ decay, the ST and DT rates no longer depend solely on the 
branching fractions of interest; the rates are modified by x, y, and the magnitudes and phases 
of various amplitude ratios. Therefore, one must correct the measured ST and DT yields in 
order to extract the branching fractions. Conversely, as shown in Ref. [96], using these yields, 
it is possible to determine simultaneously the mixing and amplitude ratio parameters as well 
as the relevant branching fractions. The estimated uncertainties on the fit parameters based on 
approximately 3x10" pairs are presented in Table 10 using efficiencies and background 
levels similar to those found at CLEO-c. 

5.6.4 Dalitz Plot Analyses 

A Dalitz plot analysis of multibody final states measures amplitudes and phases rather than 
the decay rates. A better understanding of final state interactions in exclusive weak decays is 
needed in order to elucidate the origin of CP violation in the B sector [17]. Currently, the most 
important three-body decay for 7/^3 determination is —> K^n'^n'. Recently Babar [101] 
and Belle [75] have reported 7 = (70 ± 3rJ^}|)° and ^3 = (77:}^ ± 13 ± 11)°, respectively, 
where the third error is the systematic error due to modeling of the Dalitz plot. A fit which 
includes only established resonances, similar to that of CLEO [97] which used Breit-Wigner 
(BW) line shapes, provides very poor description of the data. An improved description requires 
the inclusion of two ad-hoc tttt 5 -wave resonances which account for ~10% of the fit area. 

It may be possible to reduce the systematic error due to model uncertainty by using the 
K-matrix formalism to describe the tttt 5 -wave contribution to D'^ K'-^n'^n'. The FOCUS 
collaboration has used a hybrid BW/K-Matrix model where the higher spin resonances are de- 
scribed using BW line shapes and the nn S -wave is parameterized using the K-matrix model of 
Anisovich and Sarantsev [102] to describe D'^ —> n^n~n^ [103]. Similarly, the CLEO collabo- 
ration has searched for nn S-wave in D° n^n'n'^ [104] following the ^-matrix formalism of 
Au, Morgan, and Pennington [105]. Currently, the CLEO collaboration is considering K-matrix 
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descriptions of both the nn S-wave and the Kn P-wave contributions to — > K^^ti^tC . This 
more sophisticated description of the J^^n'Tf decay may lead to a smaller systematic error in 
the determination ofyl<p-i. 

Both D° and 5° populate the Dalitz plots K^^n^n-, n^n-n^, K^K'n^ and K^K^n'' and so 
can be used in the determination of 7/^3 which exploit the interference between b —> cus 
(B~ D^K~) and b ucs (B~ D^K~) where the former process is real and the later is 
oc e"'^ [106]. Studying CP tagged Dalitz plots allows a model independent determination of the 
relative Z)° and phase across the Dalitz plot. Consider Z)^ K'^n'^n" decay which proceeds 
through intermediate states that are CP+ eigenstates, such as K^fo, CP- such as K'^p and fla- 
vor eigenstates such as K*-n-' [97]. The Dalitz plots for tA'(3770) ^ D^D^ S+K^n^n' and 
(/'(3770) D'^D^ S ^K'^tt'^tt" will be distinct and the Dalitz plot for the untagged sample 
i/^(3770) D^^rP XK^7r'^7T~ will be different from that observed with uncorrelated D mesons 
from continuum production at ~ 10 GeV [97]. The CLEO collaboration will perform a simul- 
taneous fit to CP+, CP-, and flavor tag samples with BW/K-matrix hybrid models. A good fit 
will validate Dalitz plot model and is expected to reduce the model dependent systematic error 
on the 7/^3 measurements to a few degrees. If a good model cannot be constructed, a model 
independent result can be obtained from a binned analysis of the CP tag and flavor tag Dalitz 
plots. It is noteworthy that, although the statistical sample is smaller, the negligible nn 5 -wave 
contribution to D ^ n^n'n^ may lead to a better decay model and smaller uncertainty on 7/^3 
relative to that obtained with D K^n'^n'. 

5.6.5 Summary 

CLEO-c measurements are important inputs for each of the methods proposed to determine 
7/^3. The Gronau-London-Wyler and Atwood-Dunietz-Soni methods will benefit from im- 
proved mixing parameters and DCS parameters, respectively. The Dalitz plot method will ben- 
efit from improved models of three-body charm decay. 

The quantum correlated D^D^ system from the decay of ^(3770) provides time-integrated 
sensitivity to the mixing parameters x and y and the doubly-Cabibbo suppressed parameters Rd 
and 6. Targeted analyses will provide the first measurement of cos 6 and sensitivity to Rm and y 
competitive with 5-factory measurements. A comprehensive analysis of single tag and double 
tag yields allows simultaneous determination of hadronic and semileptonic branching fractions, 
mixing parameters, and doubly-Cabibbo suppressed parameters with sensitivity similar to the 
collection of targeted analyses. 

The decays B~ followed by a three-body decay of the D such as K^n'^n' or 

n'^n'n^ currently provide the greatest sensitivity to the CKM angle 7/^3. The precision of these 
measurements will eventually be limited by the understanding of the D K^^n^n' Dalitz plot. 
K-matrix descriptions of the nn S-wave may yield improved models of charm Dalitz plots and 
these models will be tested using the CP tagged sample of charm decays at CLEO-c. The model 
uncertainty, which is currently +10°, may be reduced to a few degrees with sufficient data. 

5.7 Be Mesons: Another Probe of CP Violation 

Contribution from R. Fleischer 

Many strategies to explore CP violation through the "conventional" charged 5„ and the 
neutral 5^,, mesons were proposed in the literature. There is, however, another species of B 
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mesons, the 5,-meson system, which consists of ~ cb and ~ be. These mesons were 
observed by the CDF collaboration through their decay B^ Jjif/i^vc, with the following mass 
and lifetime [107]: 

Mb, = (6.40 + 0.39 + 0.13) GeV, tb, = (0.46:°:j^ + 0.03) ps. (63) 

Recently, DO reported the measurement of B^ J/i/ffi^X [108], yielding 

Mb, = (5.95tJ^ ± 0.34) GeV, tb, = (0.448:;|:i^^ ± 0.121)ps. (64) 

Moreover, there is now also evidence for the decay B^ J/ifrn'^ from the CDF collaboration 
[109], allowing the extraction of 

Mb, = (6.2870 + 0.0048 + 0.001 1) GeV. (65) 

As run II of the Tevatron will provide further insights into B^ physics and a huge number 
of Be mesons will be produced at the LHC, the natural question arises whether insights into 
CP violation can also be obtained from 5,.-meson decays. In order to address this question, 
we have to consider non-leptonic B^ decays, as in the case of the exploration of CP violation 
through decays of the conventional Bu,d,s mesons. Since the B^ mesons are charged particles, i.e. 
do not exhibit mixing effects, we have obviously to rely on direct CP violation. It is well-known 
that the extraction of y from such effects suffers from hadronic uncertainties, and that an elegant 
solution of this problem is offered in the case of the 5* mesons through the "triangle approach" 
illustrated in Fig. 18 [58]. 

k-dO) 



Figure 18: The extraction of y from B- K-{D^, D'^,D1} decays. 

In the 5c-meson system, such a strategy for the determination of y is offered by the de- 
cays 5* D'^D, which are the 5c-meson counterparts of the 5* K-D modes used in the 
conventional triangle method [58], and satisfy the following amplitude relations [110]: 

^A{B^e ^ D^^Dl) = A(5+^D;D°)+A(5^^D;^D") (66) 

^^2A(B- D~Dl) = A{B'^ ^ D;D°) +A(5- ^ D-D% (67) 

with 

A(5;r ^ Z);^D°) = A(b; d;d") (68) 

A(5+ ^ D+D") = A(B; D;D'') X e^'V. (69) 

At first sight, everything is completely analogous to the 5* K'D case [58]. However, 
there is an important difference [111, 112], which becomes obvious by having a look at the 
corresponding Feynman diagrams: in the 5* D'^D system, the amplitude with the rather 



A{B+ ^ K+D^) V2A(B- 
A(B+ K+DO) = A(B- R-D^) 
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A{B+ ^ D+dO) = A{B- ^ DJD^) 



Figure 19: The extraction of y from 5* D'^iD^, D", D° } decays. 

small CKM matrix element Vut is not colour-suppressed, while the larger element Vet comes 
with a colour-suppression factor. Therefore, we obtain 

4^x^*0.4x3 = 0(1), (70) 

^ \Vcb\ a2 

and conclude that the two amplitudes are similar in size. In contrast to this favourable situation, 
in the decays 5* K'D, the matrix element Vub comes with the colour-suppression factor, 
resulting in a very stretched triangle. The extraction of y from the 5* D'^D triangles is 
illustrated in Fig. 19, which should be compared with the squashed 5* K-D triangles shown 
in Fig. 18. Another important advantage is that the interference effects arising from D", D° 
n^K~ , which were pointed out in the context of the 5* K-D case in [57], are practically 
unimportant for the measurement of BR(5^ DJ^D") and BR(5^ D*D^) since the 5c-decay 
amplitudes are of the same order of magnitude. This feature implies also that the sensitivity to 
Z)"-D° mixing, which provides a nice avenue for new physics to enter these triangle strategies 
[1 13], is considerably smaller than in the 5„ case [112]. 

Consequently, the 5* D-^D decays provide - from the theoretical point of view - the ideal 
realization of the "triangle" approach to determine y [111]. On the other hand, the practical im- 
plementation still appears to be challenging, although detailed experimental feasibility studies 
for LHCb are strongly encouraged. Using a relativistic quark model [114], which predicts the 
branching ratios for D"e^Ve, K^D^ and B^ D^D'^ in good accordance with 

experiment, the branching ratios for non-leptonic 5,. decays were recently predicted, yielding 
in fact a pattern in nice accordance with (70) [1 14]. In another study [115], also semi-leptonic 
Be decays were investigated, which give a nice testing ground. 

In addition to CP violation, there are several other interesting aspects of B^ physics. Since 
these mesons are the lowest lying bound states of two heavy quarks, b and c, the QCD dynamics 
of the B^ mesons is similar to quarkonium systems, such as bb and cc, which are approximately 
non-relativistic. However, there is an important difference: as the B^. mesons contain open 
flavour, they are stable under strong interactions. The quarkonium-like B^ mesons provide an 
important laboratory to explore exciting topics such as heavy-quark expansions (HQE), non- 
relativistic QCD (NRQCD) or factorization in a setting that is complementary to "conventional" 
weak hadron decays [116], and should provide further valuable insights into the interplay of 
strong and and weak interactions. 



A(5+ ^ D+D°) 



A(5+ ^ D+DO) 



A(B- D;D°) 



AiB- ^ D-D^) 
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5.8 Extraction of 7 at LHCb with a Combined Bs ^ Df K+ and Ba ^ D±;r+ 
Analysis 

Contribution from G. Wilkinson 

5.8.1 Introduction 

The potential of ^ D*K^, ^ D*±7r^ and B^ ^ D*;r^ decays for extracting the CKM 
angle 7 is well known [116]. An analysis based on any of these modes in isolation, however, 
suffers from the problem of discrete ambiguities, and in the case of Bd decays, of difficulties 
posed by very small interference effects. This report explains how a combined analysis of, for 
example, Bs D*K^ and Bj D-;r^ under the assumption of U-spin symmetry, circumvents 
these problems. This strategy was proposed in [1 17] and is here explored in the context of the 
expected performance in these decays at LHCb. 



5.8.2 Formalism 

From measuring the two flavour tagged decay modes, Bj ^ D^K , as a function of proper time, 
t, the CP asymmetry Acp(D^K~)(0 may be constructed: 



Acp(D:K-)(0 ^ 

This has the following dependence: 
Acp(D^K-)(0 = 



Bs ^ D+K (0 - Bs ^ D+K (0 
Bs-^D+K-(0 + B:^D+K-(0' 



Cs cos Am^t + S s sin Am j 
cosh(Ar,?/2) - Aaf, sinh(Ar,?/2)' 



where Am^ and AT, are the mass and lifetime difference between the heavy and light Bj eigen- 
states, which for the purposes of this discussion are assumed to be known. The three observables 
Cs, Ss and A^r^ can then be fitted from the data. (In doing this the full flavour untagged statis- 
tics may be used to fix Aaf,-) From performing an equivalent analysis for the D~K"^ final state 
three additional observables, Cj, S , and Aaf,, can be obtained. The observables depend on the 
underlying physics parameters in the following manner: 

CsACS) = -(+)(^)' (71) 

^ _ 2r, sin(«^, + y + (-) S,) 

1 + 

2r, cos(^, + r + (-) Ss) 



AAr.,(AAr.) = (73) 

1 + rj 

Here is the ratio of amplitudes between the interfering tree diagrams, is a possible CP 
conserving strong phase difference between the diagrams and (p^ is the CP violating weak phase 
associated with the Bj-Bs oscillations, believed to be very small in the Standard Model. It is 
assumed that (ps can be constrained from other measurements, such as those in Bg —> J/i/fcp 
decays. 
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Measurement of the six observables C^, C^, S s, S s, Aaf, and Aat^ allows r^, 5s and y to be 
determined. The same information may be extracted by making a simultaneous event-by-event 
fit to the four decay distributions. 

Exactly parallel relations hold in the system, for analysis of — > D^tt'' or Bj D**7r'^ 
decays. In this case however, the negligible width splitting between the mass eigenstates, means 
that there are effectively only four observables: Q, Q, and S4. These observables de- 
pend on rd, (f>d (measured from Bd J/i/^K*' and equal to 2/3 in the Standard Model), 6d and 7. 
Note that the value of these observables will in general be diff'erent between Bj D-n"^ and 
Bd —> D**;r'^, firstly because of the possibility of different values of and 64 for the two cases, 
and secondly because the / = 1 state of the Bd —> D**n'^ decay introduces some sign flips in 
expressions 72 and 73 (see [117] for more details). 

5.8.3 LHCb Event Yields and Performance 

LHCb has reported simulation studies of Bs DfK"^ and Bd D*-n'^ in [1 18]. The results are 
summarised here, together with estimates from a preliminary study of Bd D-n"^. 

The isolation of Bj D*K^ decays is experimentally challenging, because of the low 
branching ratio and background from the order-of-magnitude more prolific Bj D^tt"^ decay 
mode. The LHCb trigger system gives good performance for fully hadronic modes and selects 
Bs D^K"^ events with an efficiency of 30%. The n - K discrimination of the RICH system 
reduces the Bs ^t^^ contamination to ~ 10%. It is estimated that the experiment will accu- 
mulate 5.4k events per year of operation (with a year defined as 2 fb~^ of integrated luminosity), 
with a background to signal level of < 0.5. The excellent ~ 40 fs proper time precision provided 
by the silicon Vertex Locator will ensure that - provided Am, is not far in excess of expectation 
- the Bs oscillations will be well resolved, and hence the CP asymmetries can be measured. It 
is estimated that the statistical precision on y from this channel alone will be 14° for 2fb~\ 
assuming Am, = 20ps \ AFv/r^ = -0.10 and taking plausible values of y and 6s- 

The channel Bd D*-Tt^ has been investigated through the sub-decay D** D*'(^ K7r);r*, 
in which it is estimated 206k events will be accumulated per year, with a background to signal 
level < 0.3. Earlier studies using inclusive D'^ decays [119] suggest that these statistics can be 
increased still further. 

Studies are underway to investigate the feasibility of reconstructing Bd with D* 

K.-n'^n-. The preliminary indications are that 210k events will be collected each year, with a 
background to signal ratio of around 0.3. 

5.8.4 Extraction of 7 from modes in isolation 

In extracting y from the observables of a single decay mode, two problems are encountered, 
one general, and one specific to the Bd decays. 

1. The extraction of y from S s(d) and S s{d) yields 8 possible solutions. The same is true 
of calculating y from Aaf^ and Aaf,. Figure 20 illustrates this by plotting all possible 
solutions for y and 6s in the case where the true values are assumed to be y = 60° and 
6 - 60°. Therefore the study of Bd Vi-it^ or Bd XT-Jt^ in isolation results in an 
8-fold ambiguity for y. The extra observables available in the Bs system mean that in 
principle there is only a 2-fold ambiguity, but in practice measurement errors may make 
it difficult to exclude local minima coming from the additional bogus solutions. This 
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problem will be accentuated if the magnitude of Ar^ /Fs is small, hence making A^r, and 
Aaf, difficult to measure. 
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Figure 20: Valid solutions for 6^ and y for (a) the flavour tagged S s, S s observables and (b) the 
untagged Aaf, , ^af, observables. In both cases there are eight solutions (although, in the case of 
(b), these are two-fold degenerate). The solid circle indicates the true solution y = 60°, 6^ = 60°. 



2. As is clear from expressions 72 and 73, extracting y from S s(d) and 5 ^y), or A^r, and 
Aaf, require that r5(rf) be known. From comparing the CKM elements in the interfering 
diagrams it is expected that ~ 0.4 and ~ 0.02. Therefore |C.v| will be significantly 
difl'erent from 1, and measurements of and C.v will allow r., to be extracted from the 
data - indeed this is what is done in the present LHCb simulation studies. This will not 
be possible, however, for r^. Instead this parameter has to be set using external assump- 
tions [120]. These assumptions introduce a troublesome systematic error to the analysis. 

Both of these problems may be tackled by making a combined analysis of U-spin related 
and Bs modes. 

5.8.5 A Combined U-Spin Analysis of ^ D*K^ and ^ D±7r^ 

The decays Bj D^K"^ and B^ D-n"^ are identical under U-spin symmetry, ie. the exchange 
of d and s quarks. This symmetry allows the observables in both decays to be combined in 
a manner to yield certain relations, which then give y under the assumption that the a priori 
unknown hadronic contributions to the observables are identical in both channels. These un- 
knowns are the strong phases 6d and 6^, and aj and a.,, the hadronic contributions to and r^, 
defined by r^/., = Ud^sff^^, where the CKM factors, /J'™, are easily calculable. 

The following analysis follows the strategy proposed in [117]. The example plots and num- 
bers assume the scenario y = 60°, (J^^ = 60°, = I, (pd = 47° and (p^ = 0°. The experimental 
contours assume that in one year LHCb can measure 5 j and 5^ with an uncertainty of 0.014, and 
S s and 5 V with an uncertainty of 0.14 (results consistent with the performance figures quoted 
in section 5.8.3). It is also useful to assume that in the early year of operation the analysis can 
benefit from studies of Bd D-n"^ made at the 5-factories. Taking existing measurements [39], 
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and scaling to 2500 fb~^ to represent a plausible 5-factory integrated luminosity in 2008 gives 
an error on S j and Sd of 0.014. 

Using expression 72, the sine observables for the Bs and Bd channels may be combined to 
give the following exact relations: 



Us cos 6s 
cos 6a 



1 / sin(0j + y) U 5, + 5/ 
7?\sin(0, + r)/Urf + Td. 



(74) 



Us sin 6s 
ad sin 6d 



1 / cos(0d + y) 



R \cos((^^ + y) 



S s S s 



(75) 



Vd « 1 R 



Herei? = (i 

1^) ( TT7 )' "^here A is the sine of the Cabibbo angle. Because , ^ - ^ - 
and hence these relations may be exploited without the need to measure r^. In the limit of 
full U-spin symmetry, the left hand sides of equations 74 and 75 are equal to unity, and both 
relations give a determination of y. 
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Figure 21: Contours formed from (a) expression 74 and (b) expression 75 showing the 1 sigma 
contours for one year of LHCb data, and one year of LHCb data together with the measurements 
which may be available from the 5-factories by 2008. The solid circle indicates the true solu- 
tion; the inverted solid triangles indicate fake solutions from an analysis of Bg D^K"^ alone. 
Exact U-spin symmetry corresponds to a value of unity for the ordinate in both plots. 



Fig. 21 show the one sigma contours which will be obtained with one year of LHCb opera- 
tion, together with the improvement possible if 5-factory data are also included. The ambiguous 
solutions from the 'conventional' analysis are indicated. It can be seen that these bogus solu- 
tions are disfavoured by the combined analysis. There is, however, a mirror solution at - 1 80° +y 
outside the region of the plots. This possibility can be rejected either by making assumptions 
about the orientation of the unitarity triangle, or by noting that such a solution is accompanied 
by a very sizable strong phase difference, which is unlikely to be the case. 

The precision achievable on y is in general different for the two expressions. In the chosen 
scenario, relation 75 gives the best result, returning an uncertainty of cry =^^^ degrees. In the 
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one year analysis the contribution of the 5-factory data is significant. With five years of data, 
the precision improves to 5°. 

The plots also allow any biases from U-spin breaking eff'ects to be assessed. It can be seen 
that relation 75 has the steepest contour and hence exhibits the best robustness. For example, a 
20% deviation in cos Sjud cos 5d from unity leads to a 3° bias in y. 

In order to further assess the reliability the y extraction, the U-spin dependence may be 
weakened. This can be done by combining relations 74 and 75 into a single expression which 
involves either 64 and 6s or Us and a^. For example, the latter exercise yields the equation 

a, ^ ^ 1 sin2(.^^ + y) / (5 cos^(0, + y) + (5;)^sin^(0, + y) 
aa ~R sin 2(0, + y) )| (5 +)2 cos2(,^^ + y) + {S-f sin^C^^ + y) ' 

It is now possible to determine y by demanding that a, = a^, but making no assumption about 
6d and dg. With this approach, the statistical precision from five years of data is about 6 degrees. 
Again the dependence is suflaciently steep that deviations in as/a^ from unity coming from 
U-spin breaking give relatively small biases in the result. 

It should be emphazised that these analyses only make use of the flavour tagged observables 
in Bs D^K"^ and D-n^. If the magnitude of AF, is sufficiently large, then measurements 
of the untagged observables Aaf, and Aaf^ will provide additional information which will help 
further in the exclusion of ambiguous solutions, and add to the ultimate precision on y. 



5.8.6 Conclusions 

LHCb will accumulate large samples of Eg —> D*K^, Bj D**;:"^ and —> D*;r^ events. In- 
dependently each of these channels may be used to extract the CKM angle y, although in the 
case of the Bd channels this requires making assumptions about rj, the relative magnitude of the 
interfering tree diagrams. When using the flavour-tagged observables alone, this y determina- 
tion carries with it an 8-fold ambiguity, which compromises the usefulness of the measurement. 

A combined analysis of Bs —> D*K^ and Bj —> D'^n^ under the assumption of U-spin sym- 
metry allows the true solution to be isolated with only a 2-fold ambiguity. Plausible assumptions 
about the size of the strong phase diff'erence or the orientation of the unitarity triangle allow the 
remaining bogus solution to be excluded. This U-spin analysis has the further benefit of ex- 
ploiting the Bj data without the need to know r^. 

The intrinsic precision of the combined analysis is competitive with other approaches. For 
example, in the example scenario considered, a 5° uncertainty is achievable after five years of 
operation. The combined analysis does not make use of the untagged observables available 
in Bs D^K"^, which provide additional information which will improve the precision still 
further. 

The systematic error associated with the assumption of U-spin symmetry can be transpar- 
ently assessed through studying the contours associated with the measurements. In some cases 
these off'er a very robust y extraction. Furthermore, a variety of separate y determinations may 
be performed, each with diff'erent U-spin symmetry assumptions. Comparison between the 
results will help in assigning the systematic error. 

Finally, an analogous exercise can be performed from considering the U-spin related chan- 
nels Bd —> D**;r'^ and Bg —> D^^K"^. The reconstruction of the latter channel at LHCb is under 
investigation. 
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6 Summary & (near) Outlook for & y 



Contribution from A. Soni 

While the direct measurement of are now in a "matured" stage since the post-5-factory 
era of first determinations in 2001-2, past year has seen significant progress in extracting y. 
Of course, the errors are still rather large but we expect reduction of these to take place quite 
rapidly in the near future. Table 1 1 summarizes the current status and the prospects for the next 
few years for and y. 

Table 1 1 : Summary of status and prospects for /3 & y; luminosity per experiment is given in 
units of 10^55 pairs. Note current central values, jS = 0.725 & y = 69°. lEE is irreducible 
experimental error and ITE is irreducible theory error 
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0.037 
+13 + 9 + 13 
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0.028 
+9 + 6 + 7 


10 
0.02 
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0.02 

? 


0.001 
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On p the current error of 0.037[121] should decrease to about 0.02 in about 3 years but 
after that the currently estimated irreducible experimental error (lEE), which both BABAR & 
BELLE agree to be around 0.02 (see talks by Browder and Lange), will start to become the 
limiting factor. Since the irreducible theory error (ITE) on /3 is estimated to be ~ 0.01[l](and 
talk by Mannel) it is important that effort is put into reducing the lEE so that it is comparable to 
ITE. 

Regarding y, so far the B DK with Dalitz analysis of the multi-body mode — > K^n^n' 
proposed by [17] and independently in [16] has been very successful. Using about 3 x 10^ B- 
pairs BELLE gets y = 68 +15 + 13 + 11° (see talk by Abe), where the first error is statistical, 
the second is experimental systematics and the third is modelling error of the Dalitz study. 
BABAR's combined use of GLW[58, 46], ADS [57] and this multi-body mode in about 2 x 10^ 
5-pairs leads to y = 70 + 31 + 12 + 13° (see talk by Schune). Combining these we arrive at 
the current error on y in the Table 1 1 of about +13 + 9+ 13. We think it is a bit too premature 
to add the modelling and the systematic errors in quadrature to the statistical error; especially 
the modelling error ought to be left separate for now. Prospects for appreciably reducing the 
error on y in the next three years or so seems to be excellent; there are many reasons for this 
optimism. In fact is well known there are many common modes CPES(GLW), single-Cabibbo 
suppressed [67] (CPNES & CPES) as well as double Cabibbo suppressed ones [57] (CPNES 
& CPES). With increasing luminosity from about 10* to about 10^ B-pairs as more and more 
of these channels get added to the analysis, it is anticipated that the statistical error on y will 
reduce faster than 1 / V^. This is the basis for guess-estimating the error in 2006 to get around 
+9 + 6 + 7. As we enter the era of 1 ab~Vexpt around 2008 then we expect the error on y to be 
around 5 to 2 degrees. Input from charm factories can help a lot in reducing the error on y (see 
talks by Asner, Atwood and Petrov). The optimistic error of 2° may well be plausible by 2008 
if sufficient information from charm factory is forthcoming. 

Another important way to get these angles is, of course, by studying time-dependent (i.e. 
MIXCP) CP violation via £)0(*)«^" Qnce again, all the common decay modes of D° 
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and D'^ can be used just as in the case of direct CP studies involving B- decays. Therefore, 
needless to say input from charm factory also becomes desirable for MIXCP studies of 
£)0(*)"^0" jj. |g important to stress that this method give not only the combinations of the angles 
(2/3 + y = a -/3 + n) but also in addition this is another way to get jS cleanly [69, 14, 1 1]. In fact 
whether one uses 5* with DIRCP or - B^ with TDCP these methods are very clean with (as 
indicated above) the ITE of « 0.01. However, the TDCP studies for getting y (with the use of jS 
as determined from ij/Ks ) is less efficient than with the use of DIRCP involving B-. However, 
once we go to luminosities > 1 ab~^ then the two methods for y should become competitive. 
This method for getting p is significantly less efficient than from the if/Ks studies [69]. 
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